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1 =   INTRODUCTION 

The  problem  of  reflection  of  electromagnetic  waves  from 
natural  surfaces  is  of  interest  in  a  number  of  different  fields . 
In  particular,  quite  a  few  attempts  have  been  made  recently  using 
the  results  of  monostatic  and  bistatic  radar  reception  of  sur- 
faces echoes  to  map  the  surface  of  terrestrial  bodies  [28,29].  In 
these  experiments  information  about  the  surface  to  be  mapped  is 
obtained  by  comparing  the  target's  echo  signal  with  the  hypo- 
thetical signal  which  would  have  been  received  from  an  ideal  tar- 
get.  The  data  is  then  interpreted  to  determine  the  electrical 
and  statistical  properties  of  the  target  surface. 

When  an  electromagnetic  wave  is  incident  on  a  plane  inter- 
face (of  sufficiently  large  dimensions)  between  two  media,  it  is 
reflected  specularly  according  to  well  known  laws.   The  reflected 
field  depends  on  the  wavelength,  the  angle  of  incidence  and  the 
electromagnetic  properties  (permittivity,  permeability  and  con- 
ductivity) of  the  two  adjoining  media.   In  fact,  the  laws  of 
reflection  by  a  plane  boundary  are  so  well  understood  that,  con- 
versely, the  electrical  properties  of  a  material  are  often 
determined  by  measuring  its  reflection  coefficient. 

On  the  other  hand,  the  incident  field  on  a  rough  surface 
will  be  scattered  into  various  directions,  though  certain  priv- 
ileged directions  may  receive  more  energy  than  others.  Rough- 
ness itself  is  a  relative  term.   The  same  surface  may  be  rough 
for  some  wavelengths  and  smooth  for  others;  or  for  the  same  wave- 
length it  may  be  either  rough  or  smooth  for  different  angles  of 
incidence.   A  smooth  surface  is  thus  the  limiting  case  of  a 


rough  one.   The  limit  depends  on  the  wavelength  and  angle  of  inci- 
dence of  the  incident  radiation. 

The  problem  of  scattering  from  a  rough  surface  has  become  of 
special  interest  in  recent  years,  particularly  in  connection  with 
the  propagation  of  radio  waves  at  frequencies  above  30  Mc/s.   A 
large  number  of  papers  [1-2,  4-6,  14-19,  40,  41]  have  been  pub- 
lished on  the  subject  of  scattering  from  rough  surfaces,  espe- 
cially in  the  last  15  years.   Extensive  experimental  data  have 
been  accumulated  and  many  theories  have  been  developed  to  explain 
and  predict  measured  data.   None  of  these  theories  are  general 
and  rigorous  at  the  same  time.   In  order  to  arrive  at  results 
that  lend  themselves  to  reasonably  simple  numerical  calculations, 
or  to  arrive  at  any  result  at  all,  certain  simplifying  assumptions 
are  introduced  into  these  theories.   By  far,  the  largest  number 
of  rough  surface  scatter  theories  are  based  on  the  Kirchhoff 
approximation  of  the  boundary  conditions  which  are  required  to 
evaluate  the  Helmholtz  integral.   Apart  from  the  original  Kirch- 
hoff  postulation,  other  methods  of  approximating  these  boundary 
conditions  have  also  been  suggested,  e.g.  estimating  the  surface 
current  distribution  from  the  incident  magnetic  field,  or  expres- 
sing the  total  field  by  means  of  local  reflection  coefficients. 
In  principle  these  approximations  are  but  different  versions  of 
the  same  Kirchhoff  method  or  the  method  of  physical  optics. 
Another  approximation  used  in  most  theories  is  to  assume  the 
surface  to  be  perfectly  conducting.   Within  the  errors  caused  by 
various  approximations ,  the  general  features  of  the  behaviour  of 
electromagnetic  radiation  reflected  by  rough  surfaces  have  thus 
become  reasonably  well  known  by  now. 


Scope  and  Outline  of  Present  Work . 

The  simplified  approach  for  scatter  from  rough  surfaces,  des- 
cribed above,  is  discussed  further  in  Chapter  2.   However,  this 
approach  is  inadequate  for  many  cases  of  practical  interest,  be- 
cause natural  surfaces  are  neither  perfectly  conducting  nor  per- 
fectly hard,  so  that  the  region  beneath  the  surface  inevitably 
contributes  to  the  scattered  signal.   Therefore,  a  model  con- 
sisting of  a  partially  reflecting  layer  and  taking  into  account 
the  sub-surface  region  will  predict  the  behaviour  of  a  natural 
surface  much  better  than  the  usual  model  consisting  of  a  perfect- 
ly conducting  rough  surface.   The  problem  thus  requires  the 
analysis  of  the  backscatter  of  electromagnetic  waves  from  a 
rough  layer  rather  than  a  rough  surface. 

The  most  general  case  of  layer  scatter  would  have  both  the 
surfaces  (boundaries)  randomly  rough,  with  no  restriction  on  the 
mean  depth.   However,  the  solution  of  such  a  general  case  will 
be  very  complex.   To  avoid  these  complexities,  certain  simpli- 
fying assumptions  are  made.   It  is  assumed  that 

(1)  the  mean  depth  of  layer  is  much  greater  than  the 
wavelength  of  the  incident  wave  within  the  layer,  and 

(2)  the  front  surface  of  layer  is  plane  and  the  back 
surface  randomly  rough. 

The  restriction  on  the  mean  depth  of  layer  is  necessary  to 
facilitate  the  calculation  of  the  field  produced  at  the  front 
surface  by  the  wave  scattered  by  the  back  (rough)  surface,  by 


considering  the  field  to  be  in  the  far  region  or  Fraunhofer  zone 
of  diffraction.   The  second  assumption  simplifies  the  analysis 
of  backscattered  field.   A  layer  having  the  front  surface  smooth 
and  the  back  surface  rough  was  chosen — rather  than  the  opposite 
case  of  rough  front  and  smooth  back — because  it  is  comparatively 
simpler  to  analyze  the  backscattered  field  in  the  former  case. 

A  backscatter  theory  for  the  simplified  rough  layer  is  de- 
veloped in  Chapter  3.   In  spite  of  the  simplifying  assumptions, 
the  considered  rough  layer  does  have  a  few  applications,  among 
others,  in  radar  astronomy.   In  radar  astronomy,  the  backscat- 
tered signal  from  the  moon  or  other  terrestrial  bodies  is  analyzed 
in  such  a  manner  as  to  yield  fundamental  information  regarding  the 
nature  of  the  target  surface.   Most  of  the  visible  side  of  the 
lunar  surface  consists  of  the  so-called  seas  which  have  a  rela- 
tively smooth,  flat  top  surface.   These  seas  are  believed  to  be 
made  up  of  some  porous,  sandy  material  which  covers  a  rough  rock- 
like  inner  crust.   Therefore,  the  lunar  surface  may  be  taken  to 
consist  of  a  layer  with  smooth  front  and  rough  back,  for  which 
the  theory  developed  in  Chapter  3  is  applicable. 

Assumptions  and  Restrictions „ 

In  developing  the  backscatter  theory  for  the  simplified 
rough  layer,  the  Kirchhoff  method  is  used  with  the  following 
assumptions  (see  also  Beckmann,  Ref .  5,  Ch.  1  &  3) : 

(1)   The  dimensions  of  the  scattering  elements  of  the 
rough  surface  are  taken  as  either  much  smaller  or 
much  greater  than  the  wavelength  of  the  incident 
radiation. 


(2)  The  radius  of  curvature  of  the  scattering  elements  is 
taken  to  be  much  greater  than  the  wavelength  of  the 
incident  radiation. 

(3)  Shadowing  effects  are  neglected. 

(4)  Multiple  scattering  is  neglected. 

05)   The  incident  wave  is  plane  and  linearly  polarized  with 

the  E  vector  perpendicular  to  the  plane  of  incidence . 

(6)  Only  the  far  zone  field  is  calculated,  but  no  assump- 
tion is  made  as  to  its  polarization.  Any  possible 
cross-polarized  component  has  not  been  set  to  zero, 
but  simply  excluded  from  the  results . 

(7)  The  rough  back  of  the  layer  is  assumed  to  be  in  inti- 
mate contact  with  an  infinitely  thick  slab  of  a  per- 
fectly conducting  material .   The  roughness  of  the  in- 
terface is  described  by  a  random  variable  correspond- 
ing to  the  interface-level  fluctuations,  and  is  assumed 
to  have  Gaussian  statistics . 

(8)  The  beamwidth  of  the  incident  wave  is  considered  to  be 
much  smaller  than  the  mean  depth  of  the  layer. 

(9)  The  material  of  the  layer  is  assumed  to  be  lossless. 

Experimental  Work. 

In  order  to  obtain  some  direct  backscatter  measurements ,  an 
experimental  study  of  the  behaviour  of  a  rough  layer  was  under- 
taken.  The  acoustic  simulation  technique  was  used.   The  layer 
material  chosen  was  plexiglas ,  since  it  was  one  of  the  few  mater- 
ials available  which  had  a  comparatively  low  reflection  coefficient 
(0.38)  and  was  not  affected  by  water.  A  high  reflection  coefficient 


would  reduce  the  power  transmitted  into  the  layer,  thus  making  the 
resultant  scatter  at  the  rough  back  too  small  to  be  easily  detected. 
The  experimental  results  are  given  in  Chapter  4.   From  the  measured 
data,  the  following  graphs  were  obtained: 

(1)  <o0>  versus  oz/X,  the  back  interface  being  rough. 

(2)  <a0>  versus  the  angle  of  incidence,  the  back  inter- 
face being  rough. 

(3)  Same  as  (1),  but  with  the  front  interface  rough. 

(4)  Same  as  (2),  but  with  the  front  interface  rough. 

where 

<o0>  =  average  value  of  normalized  radar  cross  section 

of  layer. 
oz  =  standard  deviation  of  the  heights  of  the  layer's 
rough  surface . 
X  =  wavelength  of  the  incident  wave. 
The  frequencies  used  were:  0.72  Mc,  1.00  Mc,  1.28  Mc,  1.60  Mc, 
1.90  Mc,  2.25  Mc,  3.00  Mc ,  and  3.50  Mc,   The  angles  of  incidence 
were:  0°,  5°,  10°  and  20°. 

The  experimental  results  obtained  are  discussed  in  Chapter  4. 
Only  a  partial  attempt  was  made  to  check  the  experimental  results 
of  Chapter  4  with  the  theory  developed  in  Chapter  3,  but  it  is  to 
be  noted  that  the  experiment  violated  assumptions  (8)  and  (9)  of 
the  theory,  which  makes  a  direct  comparison  difficult. 


2.   THE  GENERAL  KIRCHHOFF  SOLUTION  FOR  SCATTERING 
OF  ELECTROMAGNETIC  WAVES  FROM  ROUGH  SURFACES. 

A  large  number  of  theories  [1-2,  4-6,  14-19,  40-41]  have  been 
developed  on  the  problem  of  scatter  of  electromagnetic  waves  from 
a  rough  surface.  In  this  chapter  a  brief  summary  of  the  solution 
given  by  Beckmann  [5]  is  presented,  the  same  notation  being  used. 
In  Chapter  3,  this  solution  will  be  extended  to  cover  backscatter 
from  a  rough  layer.  This  chapter  could  not  be  made  very  short 
because  the  results  derived  here  shall  be  used  in  Chapter  3 . 

2.1.  General  Formulation  of  the  Problem . 

As  shown  in  Fig.  1,  the  Cartesian  coordinates  x,  y,  z  with 

■+->-*■ 
origin  0  and  unit  vectors  a  ,  a  ,  a   shall  be  used.   Let  the  rough 

surface  S  be  given  by  the  function 

C  =  Ux,y)  (1) 

where  £  is  a  random  variable. 

The  mean  level  of  the  surface  is  the  plane 

z  =  0. 
All  quantities  associated  with  the  incident  field  will  be  de- 
noted by  the  subscript  1  and  those  associated  with  the  scattered 
field  by  the  subscript  2.   The  medium  in  the  space  z> C  is  assumed 
to  be  free  space.   The  incident  field  E,  is  assumed  to  be  a  har- 
monic plane  wave  of  unit  amplitude : 

El  =  eJ(ki-r-cot)  (2) 

-►  -*■ 

where  |k,  |=k«2ir/X ,  and  k,  is  the  propagation  vector,  which  will 

always  lie  in  the  xz  plane,  and  r  is  the  radius  vector: 

-+■-*•     -*■-*■ 

r  =  xax  +  yay  +  zaz   .  (3) 


F1G.I.       THE    SCATTERING    GEOMETRY. 
I"  PLANE  OF  INCIDENCE. 
S=  THE  SCATTERING  PLANE.  . 

9»  ANGLE    OF    INCIDENCE. 
©a,  ©3*  SCATTERING   ANGLE. 


FIG. 2.    DERIVATION    OF     EQUATIONS  7  ). 


In  particular,  for  points  on  the  surface  S,  we  have 

-*•  -*■  -*■  ■+ 

(r)s  ■'  xax  +  yay  +  5(x,y)az   .  (4) 

The  angle  of  incidence,  included  between  the  direction  of 

propagation  of  ll^  an(^  the   z  axis,  will  be  denoted  by  ^i  .  The 

-►      -»■ 
scattering  angle,  included  between  az  and  k2,will  be  denoted  by 

&2r   with  8^  and  62  measured  in  opposite  senses  from  the  positive 

z  axis  (Fig.  1)  <,   Here  k2  is  the  reflected  propagation  vector,  lying 

in  the  xz  plane  with 

|k2|    =    |kj    =   k   =    2ir/X    . 

-+   -*■ 
For  lateral  scattering  out  of  the  plane  of  incidence  (k^,  az) ,  a 

further  angle  63  is  introduced  (Fig.  1). 

The  time  factor  exp(-juit)  will  hence  forward  be  suppressed. 

The  polarization  of  E^  shall  be  termed  vertical  if  E^  lies 

•*■       -*■ 
in  the  plane  of  incidence  (k]_,  az)  ,  and  again  the  scattered  field 

E-  will  be  called  vertically  polarized  if  it  lies  in  the  scattering 
plane  (k2,  az) •   Similarly,  in  horizontal  polarization,  E^  and 
E2  are  normal  to  the  incidence  plane  and  the  scattering  plane 
respectively e   The  quantities  associated  with  vertical  polarization 
will  be  denoted  by  the  superscript  "+"  and  those  associated  with 
horizontal  polarization  by  the  superscript  "-" 

The  scalar  value  of  the  scattered  field  E2  at  the  point  of 
observation  P  is  given  by  the  Helmholtz  integral 


E2<p)  =  ^te-^)sas        (5' 
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where 

t  =  e*?T  (6) 

and  r'  is  the  distance  from  P  to  a  point  r  on  the  surface  S. 

No  assumption  is  made  as  to  the  polarization  of  E2(P);  any 
possible  cross-polarized  component  has  not  been  set  equal  to  zero, 
but  simply  excluded  from  the  present  range  of  interest. 

In  order  to  deal  with  plane  scattered  waves  rather  than 
spherical  ones,  we  let  R%»,  i.e.,  we  remove  P  to  the  Fraunhofer 
zone  of  diffraction;  then  as  will  be  seen  from  Fig.  2, 

k2R'  =  k2R0-k2-r  (7) 

where  Rq  is  the  distance  of  P  from  the  origin,  so  that 

f  =  eJk2Ro-Jk2-r 


(8) 


Ro 


In  Eq.  (5) ,  E  and  ill  are  the  field  and  its  normal  derivative 

3n 

on  S.   The  exact  value  of  these  two  quantities  is  in  general  un- 
known and  the  Kirchhoff  or  physical  optics  method  consists  essen- 
tially in  approximating  the  values  of  E  and  HL  on  S  and  then 

3n 

evaluating  the  integral  (5) . 

In  the  present  case  the  field  at  any  point  of  the  surface 
shall  be  approximated  by  the  field  that  would  be  present  on  the 
tangent  plane  at  that  point.   The  radius  of  curvature  of  the 
irregularities  on  the  surface  is  assumed  large  compared  with  the 
wavelength  of  the  incident  field.   Within  this  approximation  the 
field  and  its  normal  derivative  on  S  will  be: 

(E)s  =  (l+RjEi  (9) 

/3E\  .  (VE-n)s  =  j(l-R)Eiki.n  (10) 
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?&*) 


6    IS  THE   "LOCAL"  ANGLE  OF  INCIDENCE  WITH  RESPECT  TO  THE 
NORMAL  E  AT  THE  CONSIDERED   POINT. 

G,  IS  THE  "OVERALL"   ANGLE    OF    INCIDENCE  DEFINED  WITH 
RESPECT   TO  £a  ;   IT*  IS   CONSTANT    FOR    THE    WHOLE    SURFACE 


FIG.  3.       THE   "LOCAL"  SCATTERING     GEOMETRY. 
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where  n  is  the  normal  to  the  surface  at  the  considered  point  and 
R  is  the  reflection  coefficient  of  a  smooth  plane.   The  reflection 
coefficient  R  depends  not  only  on  the  angle  of  incidence  and  the 
electrical  properties  of  the  reflecting  material,  but  also  on  the 
polarization  of  the  incident  wave.   The  Fresnel  reflection  coef- 
ficient for  a  smooth  plane  are 

+   Y^os  9-/(Y2-a  sin2e)  (11) 

R  =  Y2cos 8+/(Y2-a  sin2e) 

cos0-  /(Y2-a  sin26)  ^2) 

cos 6+  ^(Y2-a  sin* 6) 

where  a  =  1/y2  and  6  is  the  local  angle  of  incidence  (Fig.  3). 
The  quantity  Y  is  the  normalized  admittance  of  the  medium  in  the 
space  z<£  and  is  given  by 

Y  »  er/yr  (13) 

where  er  is  the  relative  complex  permittivity: 

er  ■  £/e0+  j6oXa  (14) 

with  e  the  dielectric  constant,  o  the  conductivity;  and  \ir=\i/\i0 
is  the  relative  permeability  which  may  have  an  imaginary  part  if 
the  medium  has  magnetic  losses.   The  free  space  electromagnetic 
constants  are  denoted  by  e0  and  m0 . 

2.1.  Surface  Rough  in  One  Dimension . 

Consider  a  surface  rough  in  one  dimension  only.   This  surface 
will  be  constant  along  the  y  coordinate,  so  that  n  will  always  be 
in  the  plane  of  incidence  (x,z) .   Then 
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C(x,y)  =  C(x)  (15) 

and  6  m   Bi->&   ■  e^-arctan  s'(x).  (16) 

Substituting  (8)  ,  (9)  and  (10)  in  (5)  we  get 


.   jkR 


-*■■*■       -*■        TV  •  r 

(Rv-p)»  n  eJ   ds  (17) 

S 


4irRo 
where 

*  ■*   *       r  *  *  1 

v  =  ki~k2  =k  (sin0i-sine2)ax-(cose^+cos62)az 

=  vx*x  +  vz*z  (l7a) 

♦  ■*    <+■      r  ■*    .  »■►  1 

p  ■  ki+k2=k   (sinei+sine2)ax+(cos62-cos6i)az 

-»■        -*      -*■ 
n  =  -sin6ax+cos6az 

-*■-*■  •*■ 

r  ■  xax+5 (x)az 

ds  =  secBdx 

tang  =  5 (x)  . 

For  a  surface  extending  from  x=-L  to  x=+L,  (17)  reduces  to 


E2  =  ^Ro 


kejkRo   *"   (ar-b)e^(vxx+vz5)dx 


(18) 


-L 


where 

a  -  (1-R)sin61+(1+R)sine2  (19) 

b  =  (l+R)cose2-(l-R)cose1   .  (20) 

Define  a  scattering  coefficient: 

P  =  fi  (21) 

E20 
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where  E2q  is  the  field  reflected  in  the  direction  of  specular  re- 
flection (62=61)  by  a  smooth,  perfectly  conducting  plane  of  the 
same  dimensions  under  the  same  angle  of  incidence  at  the  same 
distance,  when  the  incident  wave  is  horizontally  polarized. 
For  a  smooth  perfectly  conducting  plane  (18)  reduces  to 


"20       wRq 
From  (18) ,  (21)  and  (22) 


ikeJkR°Lcosei 


(22) 


P  ■ 


4Lcos6i 


+L 


-L 


(a£'-b)e 


j (vxx+vzC) 


dx 


(23) 


The  integral  (23)  is  easily  evaluated  when  a  and  b  are  con- 
stants, but  not  in  other  cases.   One  way  of  solving  the  general 
case  is  to  average  R  over  the  surface,  making  a  and  b  independent 
of  x.   The  other,  more  important  case  is  Y-*»  or  a  perfectly  con- 
ducting surface. 

Equations  (11)  and  (12)  then  give 

R+  =  1    ,   R"  =  -1 
so  that  a  and  b  in  (19)  and  (20)  are  independent  of  x.   Integrat- 
ing by  parts,  (23)  yields 


+  ,   ■   %        «  1+cos  (61+65)   1 

p-(e1#82)  =  +sece, *  I   J'    _. 

*■     *  xcos6i+cos62   2L 


r+L  Jv4„    et(L) 

e    dx  +  -jL±     (24) 


J  -L 


where 


+      jsec  81  sin  et    ejv«r(x) 
k(cos  61  +  cos  82 ) 


+L 


-L 


(25) 


with 


6  =82  ,      8~  -  81 
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Formula  (24)  is  the  general  solution  for  a  perfectly  conduct- 
ing, one  dimensionally ,  rough  surface. 

For  L>>\,  the  second  term  of  (24)  is  negligible  compared  to 
the  first,  so  (24)  reduces  to: 


-  -  t  11 
2L 


+l  r  * 

e    dx 


-L 


(26) 


with 


F2(6i,82)  =  sec  8^ 


l+COS(8i+8o) 


1I22L    . 


(27) 


cos8i+<Sos62 

Assume  C (x)  to  be  a  random  stationary  process  with  a  mean 
value  <$(x)>   =  0,  where  the  angular  brackets  <  >  denote  the  mean 
value . 

Representing  the  complex  conjugate  quantity  by  an  asterisk, 
and  using  (26) ,  the  mean  square  scatter  coefficient  <pp*>  is  given 
by: 


<DD*>   =  _£ 


:pp 


4L2 


+L 


J  -L 


+L 


-L 


e3Vx(xi-X2)<eJvzUi-e2>>dXldX2  (28) 


where 


and 


Si  -  C(xx)  and  £2  =  S(x2) 


,Jvz(€i-€2) 


+• 


W(zi,z2)e 


jv8(5i-52) 


dzidzn 


—  J  - 


"  X2(vz'-Vz)  <29) 

is  the  two  dimensional  characteristic  function  of  the  distribution 

W(z1#z2) . 
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Introduce  a  separation  parameter,  defined  by 

t  =  x^-X2   o  (30) 

Substituting  equations  (29)  and  (30)  in  (28)  ,  and  after 
elementary  manipulations,  we  get 


<pp*> 


Fj2 
2L 


+L 


-L 


sDVxT    [x2(vz'-vz>  J   dT       • 


(31) 


Assume  zi  and  Z2  are  normally  distributed  with  mean  values 

2 

zero,  variances  az  and  correlated  by  a  correlation  coefficient 
C  (  t ) ;  then 


WU^Zj)  = 


1  _     exp 
27ro,2/(lIC7)" 


-  zl2~2Czlz2+z22 
2az2(l-C2) 


(32) 


The  characteristic  function  for  the  above  distribution  is 
given  by 


:2(vz'-vz>  ■  exP   -v22az2(l-C) 


(33) 


Choose  a  Gaussian  correlation  coefficient 


-T2/T2 


C(t)  =  e   '  *  (34) 

where  T  is  the  correlation  distance  for  which  C(x)  will  drop 
to  the  value  e~*.   It  is  to  be  noted  that  a  surface  with  a  given 
correlation  distance  T,  appears  to  be  randomly  rough  to  an  inci- 
dent beam  only  if  the  beamwidth  is  large  compared  to  T  (i.e.,  the 

2 
area  illuminated  is  large  compared  to  T  ) ;  for  in  the  contrary 

case  the  beam  will  not  cover  a  randomly  rough  surface  with  enough 

variations  to  justify  a  statistical  description,  but  will  cover 

just  one  or  two  irregularities. 
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Substituting  (34)  in  (33)  and  expanding  it  in  an  exponential 
series,  we  get 


2_    2 


,72m  y<2me-mx2/T2 


-v2  °z 


m=0 


m! 


For  briefness  let 


(35) 


/g  ■  vzoz   =   2iraz  (cos6i+cos62) 
X 

Substituting  (35)  and  (36)  in  (31) ,  gives 


(36) 


:pp 


r>  =  U   e-9 


2L 


4 

m! 


*+L  jvXT-mT2/T2 


m=0    0 -L 


dT 


(37) 


In  order  to  avoid  a  divergent  integral  it  is  necessary  to 
rearrange  (37)  by  noting  that 

<pp*>  =  <pxp*>  +  D{p}  (38) 

where 

D{p}  =  |p  -  <p> |2  (39) 

is  the  variance  of  the  complex  variable  p . 

This  variance  D{p}  may  be  written  in  terms  of  an  integral 
over  the  scattering  surface  as 


D{P>  = 


F22 
2T 


+L 


-L 


>JVXT  [x2(vz,-Vz)  -  X(VZ)X*(V2)1 


dx   (40) 


where 


<ejvz^i>  =  <eJvzC2>  =  <ejv  ^  =  <eJvzz> 
W(z)  e:,VzZdz  =  x(vz) 


(41) 
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is  the  characteristic  function  associated  with  the  one  dimensional 
probability  function  W(z)  for  the  variable  z. 

Since  z^  and  Z2  are  normally  distributed,  therefore  z  is  also 

2 
normally  distributed  with  the  same  variance  oz  and  its  distribution 

function  is  given  by 


-*2 


w(z)    =   7-7TJ7T 


zV2oz 


(42) 


The  characteristic  function  associated  with  the  above  distri- 


bution is 


x(vz)  =  exp(-ijaz2vz2)=  exp  (-*sg)  . 


Substituting  (41)  in  (26),  we  get 


where 


<P>  =  x(vz)p0 


,  sinvxL 
~v"~L~ 


(43) 


(44) 


(44a) 


Replace  the  limits  of  integration  in  (37)  and  (40)  from  +L 
to  +00;  this  is  permissible  since  the  integral  receives  signifi- 
cant contributions  only  from  the  region  near  x  =  0.   Substituting 
(35),  (36),  (43),  (44)  in  (38)  and  evaluating  we  get 


<pp*>  =  e 


-g 


p  2+/*TF2 

2L     m=1mT7m 


gin   e-vx2T  /4m 


Note  that 


*      *.     2      <^2E2>   . 


:pp 


P  RMS   JlJo" 


(45) 


(45a) 


2.2.  Surface  Rough  in  Two  Dimensions. 

Repeating  the  procedure  of  one  dimensional  rough  surface, 
the  scatter  coefficient  p  for  a  surface  £(x,y)  rough  in  two  di- 
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mensions  is  given  by 


where 


P  = 


>2      = 
J20 


4XYCOS6- 


+X 


-X 


+Y 


(aCx+cCy-b)ejV°rdxdy   (46) 


-Y 


jke^^Acosei 
E20  =  ^-J^rT 


(47) 


A  =  4XY  is  the  area  of  S  projected  into  the  xy  plane. 


.„[, 


I  (sin0i-sin62cose3)ax-sin62sine3ay-(coseifcose2)azj 
=  vx^x+vy^y+vz^z  ^48^ 

a  ■  (l-R)sin81+(l+R)sine2cose3  (49) 

b  =  (l+R)cose2-(l-R)cos6i  (50) 

c  =  (1+R)sine2sine3  (51) 

For  a  perfectly  conducting  surface 


R+  -  1   ,   R"  =  -1 

and  a,  b,  c  are  constants. 

Integrating  by  parts  (46)  reduces  to 


p-(elf.e2fe3)   -  +  -i 


+x 


f+Y 


J  -X  J  +Y 


e3V*rdxdy  +ei2LXL 


where 


F,  - 


l+cos6icos6  2-sin8isin62cos6  3 
cose^ (cose^+cos82) 


e(X,Y)  =  - 


v, 


+x 


J-x 


jv.  r 


+Y 


-Y 


dx-i* 


+Y 


-Y 


jv.r 


+X 


-X 


dy 


(52) 


(52a) 


(53) 


Equation  (53)  represents  the  edge  effect  term  and  is  negli- 


gible compared  to  the  first  term  of  (52)  when  A>>X  t   then  (52) 
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reduces  to 


,+   +  F3  I 


jv°r 


dxdy 


(54) 


Assuming  £(x,y)  to  be  a  random  stationary  process  with  mean 
value  <  £(x,y)>  =  0,  the  mean  square  scatter  coefficient  for  a 
surface  rough  in  two  dimensions  is  given  by  <pp*>  and  is  found 
by  using  (54)  : 


F32 

<pp*>  =  _± 


+X 


4   -X  J 


+X 


-X  J  -Y 


+Y 


-Y 


exp  jvx(x1-x2)+jvy(y1-y2)|  <ejv  ^1"W)>« 

dx^dX2dyidy2     (55) 


where   ^i  -  5(xi,yi)   ,  %2  ~    5(X2'Y2)« 

Introduce  a  separation  parameter  x  defined  as 

x=   [(x2-Xl)2   +  (y2-Yi)2]Js    •         (56) 
Substituting  (56)  in  (55)  and  carrying  out  elementary 


manipulations,  (55)  reduces  to 
„   *%  2irF32 


:pp 


J0(vxyT)<e         >xdx 


(57) 


vl  0 


2,^5 


where  v   =  (vx  +vv  ) 
xy     *   y 


=  k  |  (si 


2  27  Jj 

inei-sine2cos63)   +  (sine2sine3) 

(58) 


and  J0(vxvt)  is  a  Bessel  function  of  the  first  kind  and  order 
zero  with  argument  vxyx .   Since  only  the  regions  near  x  =  0  con- 
tribute to  the  integral  and  the  contribution  from  (X,Y)  to  (-,•) 
is  negligible,  the  range  of  integration  in  (57)  has,  for  conven- 
ience, been  made  infinite . 
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To  avoid  a  divergent  integral  (57)  is  substituted  in  (38) 
to  give  the  variance 


D{p}   = 


2irF3' 


Jo(vxyT)|x2(vZf-v 


*(vz)j 


-vz)    -   x(vz)x*(vz)      tdx  .  (59) 


Substituting  (35),  (36),  (41) ,  (43)  in  (59),  we  get 


D<p)  =  i!£i! 

A 


Jo 


ovvxy 


[• 


•  £  e^""2^..-, 


m=0 


in: 


xdt 


2ttf 


3_   e-g  ]>  a? 


m=l 


ten 


J0(vxyT)e"mT   /T        TdT 


(60) 


The  last  expression  was  obtained  by  the  interchange  of 
integration  and  summation. 
Since  m  >  0,  we  have 


Jo  (VxyT>e-mT2/T2TdT  =  |2  ^T2/**        . 


(61) 


Substituting  (61)  in  (60)  gives 

TrF32T2      -g  >t     gm        «"VV    2T2/4m 
D{p}  -  — * e  y  ^    iy-     e     xy 

A        m^l  m'm 


(62) 


Substitution  (62)  in  (38)  gives 


0      A    m=i  mTm  J 


Note  that 


*      *.  -   2     =  <E2E2*> 

<PP">   ~       p~„w„   '  -  -JT     ■ 


'E20|2 


(63a) 


(63b) 


and  in  (63a) , 


Po  = 


sinvxXsinVyY 


(64) 
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Expressions  have  been  derived  for  the  mean  scattered  power 
<pp*>  for  a  statistically  rough  surface  in  Eqs.  (45a)  and  (63a), 
for  roughness  in  one  and  two  dimensions,  respectively. 
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3.   BACKSCATTER  OF  ELECTROMAGNETIC  WAVES 
FROM  A  ROUGH  LAYER 

The  case  of  backscatter  of  electromagnetic  waves  from  a 
rough  layer,  with  a  plane  upper  boundary  and  a  statistically  rough 
lower  boundary,  is  considered  in  this  chapter.   The  geometry  is 
shown  in  Fig.  4.   The  numbers  1,  2  and  3  are  used  to  denote, 
respectively,  the  initial  medium  through  which  the  incident  wave 
travels,  the  layer,  and  the  medium  which  lies  on  the  other  side 
of  the  lower  boundary  (Fig.  4) .   The  terms  "front  surface"  and 
"back  surface"  shall  be  used  interchangeably  with  "upper  boundary 
and  "lower  boundary,"  respectively. 

3.0  E  Fields  in  Various  Regions 

It  is  assumed  that  mediums  1  and  2  consist  of  some  lossless 
dielectric  materials  having  the  electromagnetic  constants  e^,  ffv 
and  e2,  Vyi   respectively.   Medium  3  is  assumed  to  consist  of  a 
perfectly  conducting  material.   The  wavelength  in  medium  1  is  given 
by  A^  and  in  medium  2  by  X2. 

The  results  of  this  chapter  have  been  derived  under  the 
restrictions  of  chapter  2,  as  given  by  Beckmann  [5].   A  list  of 
the  notation  used  is  given  at  the  end  of  the  chapter. 

The  Cartesian  coordinates  x,  y,  z  with  origin  0  at  middle 
point  of  illuminated  area  on  lower  boundary  (rough)  of  layer  shall 
be  used  (Fig.  4) .   The  rough  boundary  will  be  given  by  either  of 
the  two  functions 

5  -  C(x)  (1) 

£  *  5(x,y)  (2) 
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UPPER 


BOUNDARY 


MEAN  DEPTH 
OF  LAYER "d 


„    LOWER  BOUNDARY 


So=AREA  AB 
A«  =  AREA  CD 
L  *OC*OD 


MEDIUM  1 


6a^ 

MEDIUM  2 


FIG. 4.    THE  "SCATTERING     GEOMETRY"    FOR    ROUGH    LAYER. 
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depending  on  whether  the  roughness  is  in  one  or  two  dimensions. 
The  mean  level  of  this  boundary  is  the  plane  z  =  0  and  the  mean 
depth  of  the  layer  will  be  denoted  by  d. 

It  is  assumed  that  the  source  of  the  incident  wave  is  a  radar 
situated  at  a  point  P  in.  medium  1,  such  that  the  layer  is  at  a 
large  distance  from  it,  in  the  Fraunhofer  zone  of  diffraction 
(Fig.  4) .   If  the  beamwidth  of  the  radar  is  small  enough  to  make 
sl  K<   Rl2  (Fi9*  5) '  tnen  tne  electromagnetic  wave  incident  on  the 
upper  boundary  of  the  layer  (the  plane  boundary)  may  be  considered 
to  consist  of  a  plane  harmonic  wave. 

The  incident  wave  on  the  upper  boundary  of  the  layer  is  plane 
and  linearly  polarized  with  the  E*  vector  either  in  the  plane  of 
incidence  (xz)  or  perpendicular  to  it,  and  is  given  by 

j  (k,  .r  -  ait) 
EQe  over  the  beamwidth 


Ei  = 


(3) 


0  elsewhere 

where  k±   is  the  propogation  vector  of  magnitude  k^  =  2tt/X^    (4) 

which  will  always  lie  on  the  xz  plane,  and  r  is  the  radius  vector 

r  -  xax  +  ya  +  zaz    .  (5) 

In  particular,  for  points  on  the  upper  boundary  of  the  layer, 
we  have 

r  =  xax  +  ya  +  daz    .  (6) 

So 
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A  part  of  the  energy  within  the  incident  beam  will  be  rejected 
by  the  upper  boundary  of  the  layer  and  a  part  will  be  refracted 
into  the  layer  (Fig.  4) . 

The  angle  of  incidence  included  between  the  direction  of 
propogation  of  E,  and  the  direction  of  the  z  axis  will  be  denoted 
by  9-,,  the  angle  of  reflection  by  e'  ,  and  the  angle  of  refrac- 
tion, included  between  the  z  axis  and  k3,  will  be  denoted  by  92, 
where  Jc3,  with  a  magnitude 

k3  =  2tt/A2    '  (7> 

is  the  propogation  vector  of  the  transmitted  wave  within  the 
layer. 

Angles  9^  and  92  are  related  as  follows: 


sin  9 

sin  9 


J  ■  /KS)  - n"  • 


The  beamwidth  of  the  radar  is  given  by  the  angle  a.   Then 
the  edges  of  the  beam  (half  power  points)  make  angles  9  *  a/2 
with  the  direction  of  the  z  axis. 

The  time  factor  exp(-jo)t)  will  henceforward  be  suppressed. 

The  reflection  of  the  incident  wave  from  the  upper  boundary 
will  be  in  the  specular  direction,  since  the  boundary  is  assumed 
to  be  plane.   Therefore,  this  reflected  energy  will  not  be  received 
at  the  radar,  except  in  the  case  of  normal  incidence.   On  the 

upper  boundary,  the  reflected  electric  field  E2  will  be  given  by 

1  ft  ■*" 
]k2.r 

E2  =  V12EQe  (9) 

where  V12  is  the  reflection  coefficient  at  the  upper  boundary  for 
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a  wave  travelling  from  medium  1  into  medium  2,  and  is  given  by  (11) 
or  (12) ,  chapter  2. 

■*>  ■*  ,   _ 

r  is  given  by  (5),  and  k2 ,  with  a  magnitude 

k2  =  2tt/X1     ,  (10) 

is  the  propogation  vector  of  the  field  reflected  at  the  upper 
boundary. 

The  electric  field  E3  of  the  refracted  wave,  just  below  the 
upper  boundary  is  given  by 

o+  -+■ 
cos  (jii        Jk3-r 
E3  -  cos  ♦  *  D12E0e  (ID 


where 


D12  =  1  +  V12 


is  the  transmission  coefficient  at  the  upper  boundary  between 

media  1  and  2  for  a  wave  travelling  from  medium  1  into  medium  2, 

k-,   is  defined  by    (7)  ,    and  the   term  cos   £l  =  ^L   (See  Fig.    5)    takes 
J  cos   <J>2        A'B 

into  account  the  decrease  in  intensity  of  Eg  due  to  broadening  of 

the  incident  beam  upon  refraction.   Fig.  5  defines  <fr^  and  <fr2. 

The  cross  sectional  area  of  the  beam  CD'  at  the  lower 

boundary  of  the  layer  will  be  greater  than  the  cross  sectional 

area  A'B  at  the  upper  boundary,  since  it  is  a  divergent  beam. 

Therefore,  the  intensity  of  the  electric  field  E4,  incident  on  the 

lower  boundary  of  layer  will  be  reduced  by  the  factor 

A'B   AB  cos  $2 

CD1"  =  CD  cos  d>   '  as  comPared  to  tne  intensity  of  E3.   Then,   E4 

is  given  by 
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rfF^rmfWw 


/ 


R,=  PN,       R  =  PM,       d=OQ,       a.*QN, 

?«NM,       fsOM,       R2=0N. 

S,=  AREA   CORRESPONDING    TO      AB'.  (EQ.I4) 


S2=AREA     CORRESPONDING     TO     CD.(EQ.I5) 


FIG.5.     FURTHER    DETAILS     OF     ROUGH    LAYER     SCATTER 


COS     (fri 

4     '   cos    (J>2 


AB  cos    <(>2 
CD  cos    <J>3 


D12Eoe 


jk^.r 


=  *g  cos  H-       D12Ene 
CD  cos  <t>3    12  ° 


* 


.r 


=  §£  D12Eoe 


jK,.r 


where 
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(13) 


Si  =  AB  COS  f«. 
S,  ■  CD  cos  $3. 


(14) 
(15) 


r  is  given  by  (5)  ,  and  angle  $3  is  defined  in  Fig.  5. 

The  beam  with  electric  field  E4  incident  on  the  rough 
boundary  will  be  scattered  in  all  directions  within  the  layer 
due  to  the  roughness  of  the  boundary.   The  scattering  pattern 
will  depend  mainly  upon  the  statistical  properties  of  the 
rough  boundary. 

Within  the  layer,  the  mean  backscattered  power  <E5E,-*> 
at  AB  on  the  upper  boundary,  which  is  assumed  to  be  in  the 
Fraunhofer  zone  of  diffraction  or  the  far  zone,  is  given  by 
(45b)  or  (63b)  of  chapter  2. 


<E5E5*>  =  |E50|2Prms2  =  |E50|  <PP*> 


(16) 


where  E5  is  the  backscattered  electric  field  at  AB  within  the 
layer  and  E5*  is  its  complex  conjugate. 

For  one  dimensionally  rough  surface,  £  =  5 (x) ,  the  value  of 
E50  is  given  by 
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3*3*2 

Si  D-i^E-.koL  cos  6? 

E50  "  i   Sj  "    ,R2 "  e  '«> 

and  PRmg  is  given  by  (45),  chapter  2. 

For  the  case  of  backscatter,  which  is  being  considered  here, 
the  values  of  vx,  F2 ,  g  and  pQ  in  (45),  chapter  2,  are  given  by 
substituting  the  appropriate  values  (02  =  -8,,  6^  =  02  in  that 
order;  and  X  =  A2)  in  (17a),  (27),  (36),  (44a),  chapter  2.   Then 

vx  =  i£  sin  62  (18) 

F2  =  l/cos292  (19) 

/^T  =  r— Z  COS  6,  (20) 


sin  vYL 
v„L 


Po  =  r;:  x  tiu 


R2  =  ON  (Fig.  5)  .  (22) 

L  =  OC  =  OD  (Fig.  4) .  (23) 

T  =  Correlation  distance  of  rough  surface  (24) 

(lower  boundary  of  layer) . 
For  the  two  dimensional  rough  surface,  5  =  £(x,y),  the  value 
of  Egg  is  given  by 

E50  "   '  fe  E"ga£3f   C°5    92eJk3R2  '"' 

where  AQ  is  the  area  illuminated  on  the  lower  boundary  of  layer, 

and  prms  is  given  by  (63),  chapter  2.   The  values  of  v   ,  F3,  g 

and  pQ  in  (63),  chapter  2,  are  given  by  substituting  the  appropriate 
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values  (62  ■  -8w  6,  =  62  in  that  order;  63  =  0  and  X  =  A2)  ^n 
(58),  (52a),  (36),  (64),  chapter  2.   Since  the  backscattered  power 
is  being  calculated  in  the  plane  of  incidence,  the  values  of  vxv, 
F3,  g,  pQ  and  R2  are  given  by  (18),  (19),  (20),  (21)  and  (22), 
respectively.   L  and  T  are  defined  by  (23)  and  (24)  ,  respectively. 

The  mean  backscattered  power  <E5E5*>  at  AB  within  the  layer 
is  given  by  (16) ,  and  its  square  root  will  give  the  magnitude  of 
the  electric  field  at  that  point.   It  has  been  assumed  that  AB  is 
quite  small  and  is  situated  in  the  Fraunhofer  zone  of  diffraction. 
The  backscattered  electromagnetic  wave  over  this  region  may  there- 
fore be  assumed  to  be  a  plane  wave  given  by 

-*■     ■*■ 
jk4.r 
<E5>=  |E50|  Prms  e  (26) 

where  ppj^s  is  given  by  either  (45)  or  (63),  chapter  2;  r  is 
given  by  (6) ,  and 

£4  =  -ic3  (27) 

is  the  propogation  vector  of  the  backscattered  wave  at  AB. 

The  exact  phase  of  the  backscattered  wave  at  AB  will  vary 
from  point  to  point.   In  (26)  the  phase  has  been  approximated  by 
its  mean  value  exp(jic4.r).   This  may  be  done  because  AB  is  quite 
small  compared  to  d  and  is  in  the  far  zone  of  diffraction. 

Then,  the  mean  backscattered  field  <Eg>  at  AB  on  the  upper 
boundary  of  the  layer,  in  medium  1,  will  be  given  by 

<Eg>  =  D21  |E50iPRMS  e  (28) 

where  Prms  is  given  by  either  (45a)  or  (63a),  chapter  2;  r  is  given 
by  (6); 


k5  =  -kx 
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(29) 


is  the  propogation  vector  of  the  backscattered  wave  in  medium  1, 
and 


D21  =  1  +  V21 


=  1  -  V 


12 


(30) 


is  the  transmission  coefficient  at  the  upper  boundary  for  a  wave 
travelling  from  medium  2  into  medium.  .1. 

3.1,  Radar  cross  section  of  rough  layer 

The  radar  cross  section  a  of  a  target  is  defined  as  the 
area  intercepting  that  amount  of  power  which,  when  scattered 
equally  in  all  directions,  produces  an  echo  at  the  radar  equal 
to  that  from  the  target.   Mathematically, 

Power  reflected  towards  radar  receiver/unit  solid  angle 
0  ~  incident  power  density  on  upper  boundary  of  layer/4  it 


lim.  4irR' 
r->.» 


(30a) 


where 

R  =  distance  between  radar  and  target, 

Er  =  reflected  electric  field  strength  at  radar  receiver, 

E^   =  strength  of  electric  field  incident  on  target. 

For  the  layer,  the  mean  value  of  radar  cross  section  is  given 

by 


<o>   ■   lim  4irRn 


Ri 


<E(P)  > 


(31) 
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where 

R  =  PN  (Fig.  5)  . 

<E(P)>  =  Mean  value  of  the  backscattered  electric  field  at 
the  radar  receiver. 

It  is  more  convenient  to  calculate  the  mean  value  of  aQf 
the  normalized  radar  crioss  section  (or  differential  cross  section) 
of  the  target,  which  is  defined  as 


*'k 


lim.  — =J— 
R^-   so 


E(P) 


El 


(33) 


where  S  is  the  area  of  the  illuminated  surface  AB,  on  the  upper 
boundary  of  layer. 

In  (33)  the  only  unknown  quantity  is  E(P) .   Before  E(P)  is 
evaluated,  the  incident  wave  (E^)  will  be  redefined  so  as  to  be 
in  conformity  with  actual  practice.   Instead  of  the  incident  wave 
being  considered  a  continuous  wave,  henceforward  it  shall  be 
assumed  to  consist  of  a  pulse.   If  the  pulse  is  long  enough,  then 
it  is  a  good  approximation  to  a  continuous  plane  harmonic  wave 
while  it  lasts.   Pulsed  incidence  is  of  two  types:   (1)  Beam- 
width  limited  (2)  Pulse-width  limited.   In  what  follows,  E±   shall 
be  assumed  to  consist  of  a  pulse  of  the  beam-width  limited  type. 
It  will  also  be  assumed  that  there  is  no  elongation  of  the  re- 
flected or  scattered  pulse. 

The  length  of  the  pulse  is  taken  as  TQ  seconds,  and  the 

time  taken  by  the  pulse  front  to  travel  from  AB  at  the  upper 
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boundary  of  layer  to  the  lower  boundary  and  back  again  to  the 
upper  boundary  as  2tQ  seconds.   It  is  assumed  that  t  =  0,  when 
the  pulse  front  passes  through  N,  the  middle  point  of  AB. 

To  evaluate  E(P),  the  Helmholtz  integral  defined  by  (5), 

9E 
chapter  2  will  be  used.   In  this  integral  E  and  —  are  the 

electric  field  and  its  normal  derivative  on  AB  (or  surface  SQ) 

at  the  upper  boundary  of  layer  in  medium  1.   The  value  of  (E)s 

and  (**L\       will  vary  with  time,  having  different  values  in  dif- 


Ul 


so 
ferent  time  intervals.   These  values  will  depend  upon  the  pulse 

repetition  frequency  (PRF)  and  the  value  of  TQ  compared  with  2tQ. 
It  is  assumed  that  the  PRF  chosen  is  such  that  no  two  pulses 
interfere  with  each  other,   i.e.,  the  time  interval  between  two 
pulses  is  greater  than  the  total  time  taken  for  one  pulse  to 
travel  from  the  radar  to  the  layer  and  back  again.   Two  possi- 
bilities for  the  comparative  values  of  TQ  and  2tQ  are: 

(1)  T0  <  2tQ,  and 

(2)  TQ  >  2tQ  • 

When  T  <  2tQ,  the  pulse  shall  be  denoted  "short  pulse"  and  the 
values  of  E(P)  and  oQ  calculated  for  each  time  interval  shall  be 
grouped  under  Case  I.   Similarly,  when  TQ  >  2tQ,  the  pulse  shall 
be  denoted  "long  pulse"  and  the  values  of  E(P)  and  aQ   calculated 
for  each  time  interval  shall  be  grouped  under  Case  II. 

Case  I:   Short  pulse  (T0  <  2tQ). 


For  each  pulse  the  values  of  (E)e  and  [§£}„  can  be  grouped 

bo     Vn/so 

into  three  different  time  intervals.   These  values,  will  of  course 
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be  repetitive  for  every  pulse. 
(1)  Time  interval:  0  <  t  <  TQ 

During  this  time  interval,  the  total  electric  field  on  AB 
will  be  given  by  the  sum  of  (3)  and  (9),  which  represents  the 
incident  and  reflected  waves,  respectively: 

(E)s   =  [EX  +  E2]s 


r    i&     ± 


Eoe 


+  V12Eoe 


Dk9.r 


j  s0 


(34) 


The  normal  derivative  of  E  on  AB  is  given  by 


( 


iE\ 


SnJc 


(7E    •    n)c 
&o 

jEQ        e  K-^.n  +  V12e  k2.n 

I"     j^.r  j*2.r    ] 

jE0    Le  -  V12e  J  go  k±.n 


(35) 


since 


-k1.n  =  +k2.n  =  +k5.n 


Substituting  k2  =  k5  in  (6),  chapter  2,  gives 
jk5R' 


f  =  a 

R' 


(36) 


(37) 


where  R"  =  PM  (Fig.  5) . 

From  Fig.  5,  it  will  be  seen  that 


k5R'  =  k5Rx  -  k5.r' 


where 


r'   =  r  -  R- 


Therefore 

:5 


=  kc[R,+R-,   cos    (6-1-6,)]    -  kr.r 


■5iKrK2 


and    (37)    reduces   to 


1   u2' 


36 


■ 


(38) 


(*) 


jk5 [R1+R2  cos (61-82)]    -   j£c«r 


Ri 


Js, 


(39) 


where  the  substitution  R,  =  R'  has  been  made  in  the  denominator. 
This  is  justified,  since  the  point  P  is  in  the  far  zone. 
The  normal  derivative  of  f  on  AB  is  given  by 


(Ik)s  =  <™  -s>s0-  m*MiSo 


(40) 


Substituting  (34),  (35),  (36),  (39),  (40)  in  (5),  chapter  2,  and 
simplifying  we  get 

t     < 

j(jc2-ic5)  .£ 


jk5  [r1+r2cos(  e^e^  ]    + 

E(P)  =  3EoV12e        (kl'n) 


2ttR, 


;so 


ds  (41) 


where 


k^.n  =  -kicos  8- 


2ir 


k0  ■  t—  (sin  e,  '  a„  +  cos  6/  a,) 


'2  "  \ 


1   x 


'1  °z; 


(42) 
(43) 
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kg  ■  ~  (-sin  e1  ax  +  cos  6^^  az) 


(44) 


Since  the  upper  boundary  is  smooth  6^'  =  0i  and  (43)  reduces  to 


k2  -  il  (sin  6  a*  +  cos  9-,  a,) 


'l  V 


(44a) 


Substituting  (42)  ,  (44)  ,  (44a)  in  (41)  and  setting  the 
limits  of  integration,  we  get 

jk5 [R1+R2cos (e1-e2) 1 


<E(P)  >  = 


(-j)k1E0V12  cos  ele 


2itR, 


h 


(45) 


where 


*1  = 


+xo"ao  f+Yo 


-xo-ao 


J4lTginfl  X 


iila^sine- 


eAl 


dxdy  ■  SQeAl 


-Y. 


sinifr 


(46) 


*  =  XJ  xo  sin  91 


so  =  4xoYo 


aQ  =  QN(Fig.  5)  . 


(47) 

(48) 
(49) 


Substituting  (4)  and  (46)  in  (45) ,  we  get 

j2tt  [Rj+R?cos  (6-]  -6?)  -2ar>sine-|] 
<v.(x>s>       (-j)E0v12s0coseie     xx 

<E(P)> ^-^ 


sin  j> 
(50) 


This  value  of  <E(P)>  is  substituted  in  (33)  to  give  the  mean 
value  of  the  normalized  radar  cross  section  of  the  layer  for  the 
short  pulse  case  during  the  time  interval  0<t<TQt 
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<a0>  =  4*  v122sQ2Cos2e1 


sin  $ 


* 


(51) 


In  the  above  equation  as  ^  •*  0 ,  aj-"f     +  \   f or  normal  inci- 


sing 


dence  and  fLiHi.  +  o  for  all  other  angles  of  incidence.   Therefore 

it  may  be  concluded  that,  except  in  case  of  normal  incidence,  the 
value  of  <aQ>  will  be  very  small  in  all  other  directions. 
(2)  Time  interval:  T  <t<2tQ 

In  this  time  interval,  the  electric  field  on  AB  is  zero, 
since  the  field  due  to  the  wave  scattered  by  the  lower  boundary 
of  the  layer  is  not  present: 


(E)e   =  0 


(52) 


and  therefore  during  this  time  interval 

<aQ>   =  0     .  (53) 

(3)  Time  interval:  2t0<t<2t0+TQ 

In  this  time  interval  the  electric  field  present  on  AB  is 
due  to  the  wave  scattered  by  the  lower  boundary  of  layer,  and  is 
given  by  (28) 


Dkc.r 


(E)S0=     [<V]So=     [D21'E50|PRMS      6         '      ]s 


(54) 


The  normal  derivative  is  given  by 


(ln")=  (7E-n)s  =  [JD21|E50|p 


jkc.r 


RMS 


M's, 


(55) 


Substituting  (39),  (40),  (54),  (55)  in  (5)  Chapter  2,  setting 
limits  of  integration  and  simplifying,  one  gets 
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jk5  [R1+R2cos  (63^-62)  1 
<EW>-  (-"D2llE50l°RMS   g (56) 


where 

12  " 


+xo"ao 
"xo-ao 


2TrRx  -2 


(k5.n)dxdy.  (56a) 


Substituting  (4),  (36),  (42),  (56a)  in  (56),  and  solving  the  in- 
tegral, the  result  is 

.  .»   i~  1    o    «   j=^-[R1+R,cos(e1-e2)  1 

<E(p)>  =  (-J)^llE50lPRM.c;  Soc°   !  S    *      X        '  (57) 

The  above  equation  is  substituted  in  (33) ,  to  give  the  mean  value 
of  the  normalized  radar  cross  section  of  the  layer  for  the  short 
pulse  case  during  the  time  interval  2t0<t<2tQ+T0   : 

<0  „  =  4»  D2i2|e5o12  prms2  S0cob261     .  (58) 

o       :  2P  2 

xl  Eo 

The  statistical  properties  of  the  rough  boundary  of  layer  are 
contained  in  PrmS,  and  |E50|  is  a  function  of  the  mean  depth  of 
layer  d,  the  incident  field  strength  EQ,  etc.,  and  is  given  by 
(17)  and  (25)  for  the  case  of  one  and  two  dimensional  rough  sur- 
face, respectively. 

Case  II:   Long  Pulse  (T0>2tQ). 

* 

In  this  case  also  the  values  of  (E)s  and/3E\   for  each 
pulse  can  be  grouped  into  three  different  time  intervals.   These 
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values  will  be  repetitive  for  every  pulse. 

(1)  Time  interval:  0<t<2to 

In  this  time  interval  (E)e  is  given  by  (34)  and  consequently 

Do 

<o  >,  the  mean  value  of  the  normalized  radar  cross  section  of  lay- 
o 

er  is  given  by  (51) . 

(2)  Time  interval:  2t0<t<TQ 

In  this  time  interval  the  total  electric  field  on  AB  will 
be  given  by  the  sum  of  (3),  (9)  and  (28),  which  represent  the 
sum  of  the  incident  wave,  the  reflected  wave  and  the  wave  scat- 
tered by  the  lower  boundary  of  layer,  respectively. 

{E)S0  =  lEl  +  E2  +  <E6>]S0 


jk-i.r       3k2"r  3k5«r  -i 

Eoe     +V12Eoe     +D12E0PRMS  e      Jso     *  (59) 


Since  (59)  is  a  sum  of  (34)  and  (54),  therefore,  for  this  time 
interval,  <E(P)>  will  be  given  by  the  sum  of  (50)  and  (57) 

(-j)ev19s  cose^lf  I«i+«20b«(f1-iilra«0.lBf13 

<E(P)>   =  o   12   o  l        1 sinifr 

X1R1  * 

.  jiltRi+RoCosOi-e,)] 

+    (-J)D2i|E50|pRMS      SpCosej^      e      X! 

>1R1 

jlltRl+R2cos^9l~e2^  1 
^(-jjSoCoseie      \x  r  j4i  aoSine-L 

xjjfj [E°Vl2e     X  ^f1 

+D2l'E50lpRMS  1  *  (60) 
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Substituting  (60)  in  (33)  gives  the  mean  value  of  the  nor- 
malized radar  cross  section  of  the  layer  for  the  long  pulse  case 
during  the  time  interval  2t0<t<TQt 


m   4tr   SQ  cos261 
o  T— 2 

Al 


j4n   a^sin   e-i 
»U»*1     °  1aBt*D21|EM|plaM 


2 

.(61) 


As  discussed  earlier,  sin  *  ■*  1  for  normal  incidence  and  sin  <fr  *  0 

for  all  other  angles  of  incidence.   Therefore,  except  in  case  of 
normal  incidence,  the  second  term  of  (61)  will  be  much  greater  than 
the  first  term.   Again  p^g  contains  the  statistics  of  the  rough 
boundary  and  |E50|  is  a  function  of  d,  EQ,  etc.,  given  by  (17)  and 

(25)  for  the  case  of  one  and  two  dimensional  rough  surface,  respec- 
tively. 

(3)  Time  interval:  T0<t<T0+2tQ 

In  this  time  interval  (E)  is  given  by  (28)  and  subsequently 

so 

<a   >.    the  mean  value  of  the  normalized  radar  cross  section  of 
o 

layer,  is  given  by  (58). 
3.2.  Summary  of  Results 
Case  I:   Short  Pulse  T0<2t0 

In  this  case  the  pulse  duration  is  less  than  the  time  of  a 
two-way  trip  within  the  layer. 
(1)  Time  interval:  0<t<To 

2 


o  2 

4irV17zSncos    e- 

<o0>    "  12Xl2 ' 


sin   <f> 


(62) 
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Layer  appears  only  as  a  smooth  plane, 

(2)  Time  interval:  T0<t<2tQ 

<oQ>  =  0 

No  backscattered  signal  from  layer. 

(3)  Time  interval:  2tQ<t<2t0+T0 


(63) 


<aQ>  = 


47rD212|E50l2p2RMSSocos2e1 

>,2P  2 
xl  Eo 


(64) 


Layer  appears  only  as  a  rough  surface. 
Case  Hi   Long  Pulse  (TQ>2t0) 

In  this  case  the  pulse  duration  is  larger  than  the  time  of  a 
two-way  trip  within  the  layer. 
(1)  Time  interval:  0<t<2tQ 


<a    >   =  4irV122S0cos281 


sin4> 


(65) 


Layer  appears  only  as  a  smooth  plane 
(2)  Time  interval:  2t0<t<TQ 


4irSQcos    6i 
<0o> —2 


,4ira  sin9i 
(V12e:ixl  sin»   +D2ilE5olpRMS) 


(66) 


Layer  appears  smooth  at  normal  incidence  and  rough  as  angle 
of  incidence  departs  from  normal. 
(3)  Time  interval:  TQ<t<T0+2t0 
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*»P21  lE50l  PRMS  Socos  91  (67) 

<o  >  =  J— - 

*1  Eo 

Layer  appears  only  as  a  rough  surface. 

For  both  Cases  I  &  II,  the  reflection  coefficient  V12  is 
given  in  Eq.  (11)  or  (12) ,  chapter  2  and  transmission  coefficient 
D21  is  given  in  Eq.  (30);  also  the  RMS  backscattering  coefficient 
PRMS  is  <?iven  bY  Ec3*  *45)  or  (63^  '  chapter  2. 

3.3.  Discussion  of  Results 
Case  I:   Short  Pulse  (TQ<2t0) 


Time  interval:  0<t<TQ 

In  this  time  interval,  the  backscattered  signal  at  the  radar 
receiver  is  the  signal  reflected  by  the  smooth  front  surface  of 
the  layer  only,  and  the  value  of  <aQ>,  the  mean  normalized  radar 
cross  section  of  layer,  is  given  by  (62).   In  this  equation  <a0> 
is  proportional  to  sinj> ,  where 

*  -  x^  Xosin6l 

9i  =  angle  of  incidence  at  front  surface  (Fig.  4). 
X.  =  wavelength  in  medium  1. 

X  =  Beamwidth  (in  the  direction  of  the  x-coordinate)  of  the 
o 

incident  wave  at  the  front  surface. 

In  (62),  as  M  +0,  the  quantity 
xo 

siniji  +  1  for  normal  incidence,  and 

sin<J>  ^  0  as  the  angle  of  incidence  departs  from  normal. 
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Hence  <oQ>+   a  certain  maximum  value  for  normal  incidence 
and  <o0>+   0  as  the  angle  of  incidence  departs  from  normal.   This 
should  be  expected  from  a  smooth  surface. 

As  a  check,  the  value  of  <a0>  is  evaluated  for  normal  inci- 

dence.   Substituting  sin<j>  ■  1  and  cos  81  =  1  in  (62)  ,  gives 

♦ 

<oQ>  =  4,  1*4  So  '  (68) 

For  a  perfect  conductor  V12  =  1,  and  (68)  reduces  to 

<„>  =  4*  ^2  (69) 

Al 

which  is  the  well  known  result  for  the  radar  cross  section  of  a 
flat  metal  plate  of  area  SQ  [33,36]. 
Time  interval:  2t0<t<2t0+TQ 

In  this  time  interval,  the  backscattered  signal  at  the  radar 
receiver  is  from  the  rough  back  surface  of  the  layer  only,  and 
the  value  of  <oQ>  is  given  by  (64) . 

(a)  Dependence  on  61:      In  (64),  oQ  is  proportional  to  cos2elf 
where  e^  is  the  angle  of  incidence  at  the  upper  boundary  of 
layer.   Therefore,  as  6^0,  cos2 6^1  and  <oQ>+  a  certain  maximum 
value.   Again  as  tv*"§"»  cos  261-^0  and  <aQ>-»-0;  this  should  be 
expected  because  Qi"^  is  the  case  of  grazing  incidence  (at  the 
upper  smooth  boundary) ,  when  V^+l  and  the  incident  wave  will  be 
fully  reflected  at  the  upper  boundary,  no  energy  being  transmitted 
into  the  layer. 
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(b)  Dependence  on  pR^g:   In  (64),  <oQ>   is  also  proportional  to 
p2RMS'  which  *s  a  function  of  g,  where 


/g  -  4ir  j^-  cos  82 


(69a) 


a  =  standard  deviation  of  the  heights of  the  rough  surface  (back 

surface)  of  layer. 
X2  -  wavelength  within  the  layer  (in  medium  2) . 
92  =  angle  of  refraction  at  upper  boundary  of  layer. 

To  analyse  the  dependence  of  <oQ>   on  P^g*  two  cases  will 
be  considered: 

(a)  g  <<  1,  and 

(b)  g  >>  1  . 

Since  /g  is  proportional  to  oz/X2  ,    these  two  cases  correspond  to 
a  slightly  rough  and  a  very  rough  surface. 

The  mean  square  scatter  coefficient  for  the  two  cases,  as 
given  by  Beckmann  [5,  p.  88],  is  shown  below. 

For  a  surface  rough  in  only  one  dimension 


P2RMS  =  <pp*>  =   e"9Po2  + 
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For  a  surface  rough  in  both  dimensions 
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For  g  <<  X  (nearly  smooth  surface) ,  it  can  be  seen  from  (70) 

and  (72)  that  the  specular  term  (e"gp02)  is  the  dominant  term, 

since  T  <<  L  or  T2  <<  A.   Therefore  in  (64)  ,  <oQ>  will  have  a 

significant  value  only  for  normal  incidence,  and  will  decrease  as 

the  angle  of  incidence  departs  from  normal.   For  g  >>  1  (very 

rough  surface) ,  (71)  and  (73)  show  that  the  specular  term  (e~9pQ  ) 

is  no  longer  present  and  the  dependence  on  the  scattering  angle 

i*2'  e3' Fi<? d,Ch.2)is reduced  to  that  contained  in  F2  or  F3 ,  vx  or 

vxy,  and  g.   This  dependence  of  p2rms  on  g  (i.e.,  on  the  surface 

roughness  oz) ,  shows  that  as  g  increases  or  in  other  words  the 

X 
surface  becomes  more  rough,  P2RMS  decreases  and  the  scattering 

becomes  diffuse.   Therefore  as  g  increases,  the  value  of  <oQ> 

in  (64)  will  decrease. 

Case  II:   Long  Pulse  TQ>2t0 

The  above  discussion  for  Case  I  (Short  Pulse)  also  holds 
for  Case  II  (Long  Pulse) ,  since  the  value  of  <oQ>   for  this  case 
is  given  by  (65)  ,  (66)  ,  (67) ,  and  these  expressions  are  similar 
to  the  corresponding  Case  I  expressions  (62)  and  (64) . 

It  should  be  pointed  out  that  the  results  of  this  chapter 
are  an  extension  of  Beckmann's  work  [5],   The  essential  differ- 
ence is  that  in  case  of  Beckmann,  the  backscattered  signal  from 
only  one  surface  (rough)  is  considered,  whereas  in  the  case  of  a 
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rough  layer,  the  backscattered  signal  from  two  interfaces  (one 
smooth  and  one  rough)  is  considered. 

3.4.  Notation 

The  notation  for  this  chapter  is  collected  here  for  con- 
venience . 

£  ■  function  representing  rough  (lower)  boundary  of  layer. 
SQ  ■  area  AB  =  area  illuminated  on  upper  boundary  of  layer. 

A  =  area  CD  =  area  illuminated  on  lower  boundary  of  layer, 
o 

L  =  OC  ■  OD  (Fig.  5) . 

aQ  =  QN  (Fig.  5)  . 

r*  =  vector  NM  (Fig.  5) . 

?  =  radius  vector  (Fig.  5) . 

$2  ■  vector  ON  (Fig.  5) . 

d  =  mean  depth  of  layer. 

r  =  pn  ■  distance  from  radar  of  middle  point  of  area  illuminated 
on  upper  boundary  of  layer  (Fig.  5) . 

r  =  ON  ■  distance  from  origin  of  middle  point  of  area  illumin- 
ated on  upper  boundary  of  layer  (Fig.  5) . 

r'  =  pm  ■  distance  from  radar  of  any  general  point  on  the  area 
illuminated  on  upper  boundary  of  layer  (Fig.  5) . 

$1  =  angle  B'AB    (Fig.    5). 

<J>2   ■  angle  ABA'     (Fig.    5)  . 

<|>3  =  angle  D'CD    (Fig.    5). 

8±   =  angle  of  incidence  at  upper  boundary  (Fig.  4) . 

e^*  angle  of  reflection  at  upper  boundary  (Fig.  4). 

62  -   angle  of  refraction  at  upper  boundary  (Fig.  4) . 
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ei, Mi  =  electromagnetic  constants  of  medium  1  (Fig.  4). 

e2fM2  =  electromagnetic  constants  of  medium  2  (Fig.  4). 

a-  =  co  s  conductivity  of  medium  3  (Fig.  4)  . 

V,0   -  the  reflection  coefficient  at  the  upper  boundary  between 

media  1  and  2  for  a  wave  travelling  from  medium  1  into 

medium  2. 
D12  =  the  transmission  coefficient  at  the  upper  boundary  between 

media  1  and  2  for  a  wave  travelling  from  medium  1  into 

medium  2. 
v21  =  the  reflection  coefficient  at  the  upper  boundary  between 

media  1  and  2  for  a  wave  travelling  from  medium  2  into 

medium  1. 
D   =  the  transmission  coefficient  at  the  upper  boundary  between 

media  1  and  2  for  a  wave  travelling  from  medium  2  into 

medium  1. 
X,   =  wavelength  in  medium  1. 
X2  =  wavelength  in  medium  2. 

E,   =  incident  electric  field  on  upper  boundary  of  layer. 
EQ  =  I e^  J  =  magnitude  of  electric  field  incident  on  upper 

boundary  of  layer. 
E2  =  reflected  electric  field  on  upper  boundary  of  layer. 
E3  =  refracted  electric  field,  just  below  the  upper  boundary 

of  layer. 
E4  =  electric  field  incident  on  lower  boundary  (rough  surface) 

of  layer. 
E5  =  backscattered  electric  field  at  upper  boundary,  within  the 

layer. 
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E5*  ■  complex  conjugate  of  Eg. 

E50  =  defined  by  equations  (17)  and  (25),  Chapter  3. 

E--  =  backscattered  electric  field  at  upper  boundary,  in 

D 

medium  1. 
<E(P)>  =  mean  value  of  backscattered  electric  field  at  radar 

receiver, 
ic,  =  propagation  vector  of  wave  incident  on  upper  boundary  of  layer 
ic2  ■  propagation  vector  of  wave  reflected  by  upper  boundary  of  layer 
k3  =  propagation  vector  of  wave  incident  on  lower  boundary  of  layer 
k\  =  propagation  vector  of  the  backscattered  wave,  within  the  layer 
jcr  ■  propagation  vector  of  the  backscattered  wave,  in  medium  1. 
PRMS  "  mean  scJuare  backscatter  coefficient  of  lower  boundary  of 

layer,  given  by  equations  (45)  and  (63),  Chapter  2. 
o_  -  the  standard  deviation  of  the  heights  of  the  rough  boundary 

of  layer. 
T  ■  correlation  distance  of  rough  boundary  of  layer. 
a     =  radar  cross  section  of  layer. 
aQ   =  normalized  radar  cross  section  of  layer  =  a/area  of  layer 

surface. 
<oQ>  =  average  value  of  aQ* 
TQ  ■  incident  pulse  length. 
2t  =  total  time  taken  by  the  pulse  front  to  travel  from  upper 

boundary  of  layer  to  lower  boundary  and  back  again  to  upper 

boundary. 
Si  -   cross-section  of  beam  incident  on  upper  boundary  of  layer. 
S2  =  cross-section  of  beam  incident  on  lower  boundary  of  layer. 


50 


4o   RADAR  CROSS  SECTION  MEASUREMENTS  FOR  A  ROUGH 
LAYER  BY  ACOUSTIC  SIMULATION  TECHNIQUES. 

An  experimental  study  of  the  behaviour  of  a  rough  layer  was 
undertaken  to  obtain  some  direct  measurements  of  its  average 
radar  cross  section.  Acoustic  simulation  techniques  were  used 
for  this  study. 

4o0  Basic  Concepts. 

The  basis  for  acoustic  simulation  is  the  well  known  analogy 
between  acoustic  and  electromagnetic  waves.   By  considering  the 
pressure  or  particle  velocity  of  acoustic  waves  as  the  analog  of 
the  electric  field  intensity  in  the  electromagnetic  wave,  one 
can  obtain  analogous  expressions  for  parameters  of  interest  in  the 
two  cases  (e.g.,  impedance,  reflection  coefficient,  velocity  of 
propagation,  etc.)   Although  polarization  effects  cannot  be  ac- 
coustically  simulated  due  to  the  fact  that  acoustic  propagation 
is  a  scalar  phenomenon  whereas  electromagnetic  waves  are  in  general 
described  by  vectors,  the  phenomena  of  propagation,  reflection, 
refraction  and  scatter  that  are  of  interest  in  radar  can  be  stud- 
ied in  the  laboratory  by  means  of  an  acoustic  simulator.  This 
simulation  is  made  more  convenient  and  economical  by  the  proper 
choice  of  frequency  and  medium  of  propagation;  ultrasonic  waves 
in  water  are  used.   It  can  easily  be  shown  that  quantitative  in- 
formation (rather  than  merely  qualitative  analogy)  can  be  gained 
by  suitable  scaling  of  frequency,  range,  impedance,  or  other  pa- 
rameters . 
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A  block  diagram  of  the  equipment  which  was  used  for  the 
acoustic  simulation  of  a  rough  layer  is  shown  in  Fig.  6,   The 
acoustic  waves  are  generated  in  a  large  water  filled  tank  by  means 
of  an  electro-acoustic  transducer,  which  is  driven  by  an  electronr 
ic  oscillator.   The  transducer  is  capable  of  being  moved  in  three 
orthogonal  directions  relative  to  the  target;  it  can  also  be  ro- 
tated about  horizontal  and  vertical  axes  to  give  any  desired  angle 
of  incidence.  A  brief  description  of  the  experimental  equipment 
is  given  below.   Details  concerning  the  calibration  and  use  of  the 
apparatus  can  be  found  in  the  KSU  Technical  Report  EE-TR-1  [39] . 

4.1o  Acoustic  Simulator  Components . 

(A)  Pulsed  Oscillator;   The  high  powered  pulse  oscillator  (PG 
650-c,  Model  2,  Arenberg  Ultrasonic  Lab.,  Jamaica  Plain,  Massa- 
chusetts) is  a  variable  frequency,  pulse  modulated  radio  frequency 
oscillator  capable  of  delivering  300  volts  peak  to  peak  into  a 

93  ohm  load  resistor.   The  oscillator  can  deliver  output  pulses 
that  vary  in  width  from  2 usee,  to  lOOysec.  over  a  frequency 
range  of  0.5  Mc/sec.  to  5.0  Mc/sec. 

(B)  Transducers;  All  acoustic  measurements  were  made  by  means 
of   piezoelectric  barium  titnate  transducers  (piston  shaped) , 
manufactured  by  Branson  Instruments  Manufacturing  ,  Stanford, 
Connecticut.   Four  pairs  of  transducers  rated  at  the  following 
frequencies  were  used:  1.0  mc,  1.6  Mc,  2.25  Mc,  and  3.5  Mc.   Each 
pair  of  transducers  consists  of  a  transmitter  and  a  receiver 
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designated  by  ZT  and  ZI,  respectively .  ,  Table  I  gives  the  diam- 
eter and  beamwidth  of  the  respective  transducers. 

Beamwidth  of  transducers. 


Frequency  in  Mc. 

Diameter  of  active 
element  in  inches. 

Holder  diameter 
in  inches. 

Beamwidth 
in  deqrees. 

1.00 

7/8 

V4 

3.1 

1.60 

1/2 

3/4 

3.8 

2.25 

1/2 

3/4 

2.8 

3.50 

1/2 

3/4 

1.9 

(C)  Transducer  Compensator;   Functionally,  the  transducers 
operate  best  at  mechanical  and  electrical  resonance.   Electri- 
cally the  transducer  appears  to  be  a  capacitor  shunted  by  a  small 
conductance .   Electrical  resonance  is,  therefore,  achieved  by 
adding  the  proper  inductance  in  parallel  to  cancel  the  total 
capacitive  reactance.   Without  such  compensation,  the  large 
capacitance  of  the  transducer  cable  will  "pull"  the  oscillator 
frequency  out  of  the  range  of  frequencies  marked  on  the  oscil- 
lator coils .  A  transmitting  transducer  compensator  obtains  this 
resonance,  so  that  the  ultrasonic  oscillator  can  see  at  any  fre- 
quency a  relatively  non-reactive  load. 

(D)  Receivers  and  Detectors:   In  Fig.  6  it  is  shown  that  the  re- 
flected acoustic  waves  are  intercepted  by  the  receiving  trans- 
ducer (which  is  positioned  next  to  the  transmitting  transducer) 
and  converted  there  into  an  electrical  pulse  which  is  then  trans- 
ferred through  an  input  attenuator,  an  input  amplifier,  a  band 
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pass  filter,  a  mid-amplifier  attenuator,  and  two  wide.-band  ampli- 
fiers to  the  block  marked  "electronics".   The  "electronics"  is 
composed  of  the  following: 

(1)  Recycle  gate  generator  . 

(2)  Amplifier  gate  generator. 

(3)  Gated  amplifier. 

(4)  Detector. 

(5)  Video  amplifier. 

(6)  Boxcar  circuit  . 

(E)  Recording  Device;   The  acoustic  return  data  was  recorded  by 
the  use  of  General  Radio  Type  1521-A  graphic  level  recorder. 
This  recorder  has  an  input  resistance  of  1000  ohms  and  is  driven 
by  a  d-c  analog  voltage  from  the  boxcar  detector.   The  voltage 
normally  ranges  from  0  to  0.8  volts,  providing  a  maximum  stylus 
deflection  of  four  inches  on  the  chart. 

(F)  Motor  Control:  All  control  of  the  scanning  mechanism  is 
normally  done  from  a  remote  motor  control  box. 

(G)  The  tank  is  6  feet  long,  5  feet  wide,  and  5  feet  deep.   It 
is  constructed  of  1/8  inch  galvanized  steel  sheets. 

(H)   Scanning  Mechanism:   The  tank  scanning  mechanism  is  built  on 
a  steel  slide  which  is  supported  by  wheels  mounted  on  two  paral- 
lel angle  iron  rails  at  the  top  of  the  tank.   Thus  the  tank 
scanning  mechanism  can  be  placed  at  any  desired  distance  from  the 
target  up  to  a  maximum  of  44  inches. 
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The  scanning  carriage  is  mounted  upon  the  steel  slide  so  that 
east-west  motion  is  given  to  the  scanning  carriage  by  an  a-c  split 
phase  motor  which  is  fixed  to  one  end  of  the  slide.  Vertical  mo- 
tion of  the  transducers  is  obtained  through  another  a-c  split 
phase  motor  supported  by  the  carriage  and  driving  a  vertical  screw 
shaft.  On  this  shaft  ride  the  transducer  mount  and  the  transducer 
vertical  scan  assumbly. 

4.2o  Target  Description. 

The  layer  target  was  made  from  a  plexiglas  sheet  48  inches 
long,  6  inches  wide,  and  2  inches  thick.   One  face  of  the  plexi- 
glas sheet  was  machined  to  give  a  one  dimensionally  rough  sur- 
face, leaving  the  other  surface  smooth,  as  shown  in  Fig.  7.   The 
purpose  of  having  the  roughness  in  only  one  dimension  (z  coordi- 
nate) rather  than  the  more  realistic  two  dimensional  roughness 
(in  y  and  z  coordinates) ,  is  to  have  a  better  control  over  the 
statistics  of  the  roughness. 

The  heights  of  the  rough  surface  were  chosen  from  a  normal 
distribution  having  a  standard  deviation  oz  =  0.05  inches.  These 
heights  were  then  arranged  (by  trial  and  error)  to  give  a  random 
curve  having  a  Gaussian  autocorrelation  function  with  a  correla- 
tion distance  T  =  0.15  inches.   This  particular  correlation  func- 
tion was  obtained  by  arranging  the  heights,  such  that  there  were 
not  many  sudden  variations  in  the  slopes  of  the  curve  thus 
generated.   The  random  curve  was  then  traced  on  the  edge  of  the 
plexiglas,  which  was  cut  by  a  shaping  machine.   Finally,  the 
machined  rough  surface  was  sampled  by  a  profilometer  (apparatus 
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for  determining  the  profile  of  a  rough  surface)  and  the  statis- 
tics redetermined  from  the  measured  heights.   The  measured  standard 
deviation  cz  and  correlation  distance  T  were  0.046  inches,  and  0.150 
inches  respectively.   The  measured  autocorrelation  function  R(t) 
(Normalized)  is  shown  in  Pig.  8. 

4.3.  Procedure  for  the  Experiment . 

The  target  was  mounted  on  supports  and  placed  flat  on  the  surface  of 
the  water  in  the  tank,  such  that  one  surface  (front)  was  within  the 
water  and  the  other  surface  (back)  was  exposed  to  air,  as  shown 
in  Fig.  6.   The  purpose  of  this  particular  configuration  was  to 
simulate  the  perfectly  conducting  medium  behind  the  rough  surface, 
since  the  reflection  coefficient  at  the  plexiglas-air  interface 
approaches  unity. 

The  backscattered  signal  was  recorded  for  two  different  con- 
figurations of  the  layer.   In  one  configuration,  which  shall  be 
termed  the  "rough  back",  the  front  surface  of  the  layer  was  smooth 
and  the  back  surface  rough.   In  the  other  configuration,  which 
shall  be  termed  "rough  front",  the  front  surface  of  the  layer  was 
rough  and  the  back  surface  smooth.   In  each  case  the  backscattered 
signal  from  the  rough  interface  only  was  recorded,  the  signal  from 
the  smooth  interface  being  gated  out. 

Transducer  Positioning;   The  transducers  were  so  positioned  that 
they  were  focussed  on  the  layer  when  placed  at  a  distance  of  32 
inches  from  it.   The  transducer  assembly  was  moved  from  one  end 
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of  layer  to  the  other,  so  that  the  full  length  of  layer  was  scan- 
ned.  The  backscattered  signal  was  recorded  by  the  graphic  level 
recorder  during  the  scan. 

Frequency  and  Angle  of  Incidence;   In  the  "rough  back"  configura- 
tion of  the  layer,  the  frequencies  used  were:  0.72  Mc,  1.0  Mc, 
1.28  Mc,  1.6  Mc,  1.9  Mc,  2.25  Mc,  3.0  Mc  and  3.5  Mc.   For  the 
"rough  front"  configuration  the  frequencies  used  were:  1.0  Mc, 
1.6  Mc,  2.25  Mc  and  3.5  Mc. 

Four  different  angles  of  incidence,  0  ,  5  ,  10  ,  and  20   (at 
front  surface  of  layer)  were  used  for  each  configuration  of  the 
layer.   The  distance  of  the  transducers  from  the  front  interface 
of  the  layer  was  always  kept  constant  at  32  inches,  as  stated 
earlier. 

Pulse  Length:  The  pulse  length  chosen  was  20ysec,  which  is  less 
than  the  time  taken  (about  30ysec)  for  the  pulse  front  to  travel 
from  the  front  surface  of  the  layer  to  the  back  surface,  and 
back  again  to  the  front  surface.  This  way  two  separate  and  dis- 
tinct signals  from  the  two  interfaces  of  the  layer  were  obtained 
at  the  receiver.  These  two  backscattered  signals  were  separated 
in  time  by  about  10  u sec  for  normal  incidence. 

The  pulse  length  was  deliberately  chosen  to  be  less  than 
30Msec,  so  that  the  backscattered  signal  from  the  rough  inter- 
face of  the  layer  could  be  recorded  separate  from  the  backscat- 
tered signal  by  the  smooth  interface. 

Recording  and  Measurements :   The  value  of  the  backscattered  signal 
at  32  inches  from  a  water-air  interface  was  recorded  for  each 
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frequency o   These  values  were  then  used  to  normalize  the  backscat- 
tered  signal  from  the  layer.   The  value  of  the  recorded  signal 
directly  gives  So~^   ,  the  square  root  of  the  normalized  radar  cross 
section  of  the  layer.   Therefore,  the  square  of  the  recorded  sig- 
nal gives  o0.   The  mean  value  and  variance  of  o0  were  calculated 
for  each  set  of  recorded  data.   Further  details  and  sample  calc- 
ulations for  o0  are  given  in  Appendix  II. 

Limitations  for  20°;   For  the  20°  angle  of  incidence,  the  back- 
scattered  signal  from  only  four-fifths  length  of  the  layer  could 
be  recorded  due  to  the  limitations- in  the  transverse  motion  of 
the  scanning  mechanism.   The  measured  autocorrelation  function 
R(t)  (normalized)  for  this  part  of  the  rough  surface  of  layer  is 
shown  in  Fig.  9:  it  is  very  similar  to  the  measured  autocorrela- 
tion for  the  full  layer  shown  in  Fig.  8.   The  mean  value  and 
variance  of  the  heights  of  the  rough  surface  are  almost  the  same 
for  the  two  cases . 

It  may  be  pointed  out  that  the  sound  wave  transmitted  into 
the  plexiglas  layer  at  an  angle  of  incidence  other  than  normal 
has  two  modes  of  propagation:  the  longitudinal  mode  and  the  trans- 
verse (shear)  mode.   At  large  angles  of  incidence  a  considerable 
amount  of  sound  energy  in  plexiglas  will  be  propagated  in  the 
transverse  mode.   Therefore,  the  measured  backscattered  data  for 
the  "rough  back"  configuration  of  the  layer  at  20°  angle  of  inci- 
dence may  not  be  very  accurate. 

4.4.  Results; 

The  results  of  the  "rough  back"  configuration  of  the  layer 

[Text  continued  on  page  66] 
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have  been  summarized  in  Pigs.  10  and  11.   Of  these,  Fig.  10  gives 
the  mean  value  of  o0  (normalized  radar  cross  section)  as  a  func- 
tion of  oz/X  showing  four  curves  for  angles  of  incidence  of  0°, 
5°,  10°  and  20°,  respectively.   Here  oz/X  is  the  ratio  of  az, 
the  standard  deviation  of  the  heights  of  rough  surface  of  layer, 
to  X ,  the  wavelength  of  the  incident  wave  in  plexiglas .   The 
various  values  of  az/X  are  given  in  Table  III,  Appendix  I.   Fig. 
11  shows  the  dependence  of  <a0>  on  angle  of  incidence,  with  fre- 
quency as  a  parameter.   The  measured  mean  values  and  standard 
deviation  of  o0  are  tabulated  in  Table  IV,  Appendix  III. 

Similarly,  Figs.  12  and  13,  summarize  the  results  for  the 
"rough  front"  configuration  of  the  layer.   In  this  case  oz/X  is 
the  ratio  of  the  standard  deviation  of  the  heights  of  rough  sur- 
face of  layer  to  the  wavelength  of  the  incident  wave  in  water. 
The  values  of  oz/X  used  are  given  in  Table  II,  Appendix  I.  Again 
the  measured  mean  values  and  standard  deviation  of  o0  are  tabu- 
lated in  Table  V,  Appendix  III. 

4.5.  Discussion  of  Results. 

The  dependence  of  <a0>  on  oz/X  and  the  angle  of  incidence  as 
shown  in  Figs.  10  to  13  seems  to  be  in  accordance  with  the  the- 
oretical results  of  Chapter  3,  as  seen  below. 

(A)   Rough  Back  Configuration  of  Layer . 

(1)  Dependence  on  oz/X:  Figure  10  shows  that  for  a  particular 
angle  of  incidence,  as  oz/X  increases  (or  in  other  words  as  the 
surface  becomes  more  rough),  the  value  of  <o0>  decreases.   This 
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should  be  expected  from  the  theoretical  results  (Eq.  64,  Chapter 
3)as  discussed  in  Chapter  3.  It  is  to  be  noted  that  in  Fig.  10  the 
value  of  <a0>  is  maximum  for  6  =  0°,  falls  for  6=5°,  but  rises 
for  6  =  10°,  and  finally  falls  again  for  6  =  20°.   In  short, 
almost  all  the  curves  in  Fig.  11  are  peaked  at  0  =  10  .  The 
probable  reason  for  this  behaviour  is  the  fact  that  the  rough 
surface  has  more  facets  with  10°  slope  than  facets  with  5°  slope. 
A  check  of  the  number  of  facets  having  5°  and  10°  slopes,  respec- 
tively, confirmed  this  fact.   The  slopes  were  calculated  by 
sampling  the  rough  surface  at  intervals  of  0.05  inches,  finding 
the  difference  in  heights  of  the  adjacent  sampled  points  and  then 
dividing  this  difference  by  0.05  inches  (the  sampling  distance). 
From  these  calculated  slopes,  the  total  number  of  5°  and  10° 
slopes,  respectively,  were  counted  and  finally  compared. 

The  theoretical  results  of  Chapter  3  indicate  that  as  the 
angle  of  incidence  6  increases,  the  value  of  <o0>  should  decrease. 
However,  these  theoretical  results  were  obtained  for  a  surface 
having  a  normal  distribution  of  slopes,  while  the  rough  sur- 
face of  the  layer  used  for  the  experimental  work  does  not  have 
a  normal  distribution  of  the  slopes,  as  is  evident  from  the  fact 
that  the  10°  slopes  are  larger  in  number  than  the  5°  slopes. 
Therefore,  there  is  reason  to  expect  some  differences  between 

the  experimental  results  and  the  theory. 

* 

(B)   Rough  Front  Configuration  of  Layer. 

(1)   Dependence  on  oz/X:   in  Fig.  12,  it  is  noted  that  in  agree- 
ment with  the  theory,  the  value  of  <a0>   decreases  with  increase 
of  oz/A  (i.e.,  as  the  surface  becomes  more  rough)  as  az/\   changes 
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from  0.787  to  lo250;  however,  for  oz/A  =  1.770  the  value  of  <o0> 
increases,  although  at  oz/A  =  2.760,  <o0>  decreases  again. 

This  peaking  of  the  curves  at  az/A  =  1.770  is  probably  due 
to  some  periodicity  in  the  rough  surface  of  layer.   As  shown  in 
Appendix  IV,  the  equation  from  grating  theory  indicates  that  if 
A/A,  the  ratio  of  the  wavelength  of  the  incident  radiation  to 
the  wavelength  of  periodicity  of  the  surface,  is  small ,  then  the 
incident  wave  will  be  broken  up  into  many  scattered  waves  or 
modes.   The  autocorrelation  function  of  the  rough  surface  of  layer, 
given  in  Fig.  8,  does  indicate  that  a  periodic  surface  has  been 
superimposed  on  the  randomly  rough  surface.   Figure  8  shows  that 
the  wavelength  A  of  the  periodic  surface  is  of  the  order  of  17.3mm. 
At  2.25  Mc  (az/A  =  1.770),  the  wavelength  of  sound  in  water  is 
A  =  0.666mm.   The  ratio  j  =  —yr"  =  j6<<lr  thus  satisfying  the 
condition  stated  above.   The  scattering  pattern  should  therefore 
consist  of  the  superposition  of  the  multi-lobed  modes  of  the 
periodic  component  of  the  rough  surface  and  the  diffuse  scatter- 
ing pattern  of  the  random  component  of  the  rough  surface.   The 
resultant  scattering  pattern  will  therefore  have  a  number  of  di- 
rectional maximas.   It  is  possible  that  the  measured  backscat- 
tered  data  at  a  /A  =  1.770  was  on  one  of  these  maximas,  whereas 
the  measured  backscattered  data  for  the  adjacent  value  of  oz/A  = 
1.25  was  not  on  one  of  these  maximas. 

Therefore,  one  probable  reason  for  the  peaking  of  the  curves 
at  az/A  =  1.770  (frequency  =  2.25  Mc)  is  the  grating  effect  of 
the  periodic  component  of  the  rough  surface. 
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(2)   Dependence  on  Angle  of  Incidence  6;   Figure  13  shows  that 
the  values  of  <a0>  at  6  =  5°  and  9  =  10°  are  comparable  for  all 
cases  except  for  the  case  of  1.6  Mc.   This  should  be  expected, 
because  the  rough  surface  has  a  larger  number  of  facets  having 
10°  slope  compared  to  facets  having  5°  slope,  as  mentioned  ear- 
lier.  In  case  of  the  1.6  Mc  curve,  there  is  a  big  increase  in 
<o0>  as  9  is  increased  from  5°  to  10°;  this  difference  in  the  be- 
haviour is  probably  due  to  the  grating  effect  of  the  periodic 
component  of  the  rough  surface  as  already  explained.  Otherwise, 
all  the  four  curves  do  have  a  tendency  of  decreasing  with  the 
increase  in  9 . 

Theory  indicates  that  the  value  of  <o0>  should  fall  stead- 
ily with  the  increase  in  6,  for  all  frequencies.   In  case  of  the 
curves  of  Fig.  13,  this  holds  except  for  2.25  Mc,  the  difference 
being  due  to  the  reasons  given  above. 
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5 .   CONCLUSIONS 

The  field  scattered  from  a  rough  surface  is  known  to  be  a 
complicated  function  of  the  parameters  involved,  e.g.  angles  of 
incidence  and  scatter,  polarization,  frequency,  dielectric  and 
statistical  properties  of  the  scatterer,  etc.   The  problem  is 
more  complicated,  however,  when  the  target  consists  of  a  rough 
layer  covering  a  thick  core,  this  being  the  simplest  form  of  the 
general  multi-layer  problem.   To  better  understand  the  scatter- 
ing behaviour  of  a  rough  layer,  both  theoretical  as  well  as 
experimental  work  was  attempted. 

In  Chapter  3  expressions  were  derived  for  the  mean  radar 
cross  section  (normalized)  <o0>  of  a  rough  layer.    The  incident 
wave  was  considered  to  be  a  pulse  of  the  beamwidth  limited  type, 
so  as  to  conform  with  actual  practice.   To  simplify  the  analysis 
the  rough  layer  was  assumed  to  have  a  smooth  front  surface  and  a 
rough  back  surface.   Two  important  restrictions  were: 

(1)  The  beamwidth  of  the  incident  wave  is  much  smaller  than 
the  mean  depth  of  the  layer. 

(2)  The  mean  depth  of  layer  is  much  greater  than  the  wave- 
length of  the  incident  wave  within  the  layer;  in  other 
words,  within  the  layer,  the  backscattered  electric 
field  at  the  smooth  front  interface  is  in  the  far  zone 
of  diff ration. 

From  the  theoretical  expressions  obtained  for  <oQ>   in  Chapter  3, 
it  is  reasonable  to  conclude  that  the  value  of  <o0>  decreases  as : 

(1)  the  angle  of  incidence  departs  from  normal, 

(2)  the  roughness  of  the  rough  surface  of  layer  increases. 
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The  limitations  of  the  particular  model  of  rough  layer  chosen 
for  the  theoretical  analysis  have  been  noted,  with  a  view  to  help 
improve  future  work.   This  indeed  is  a  first  step  towards  the  solu- 
tion of  a  more  general  rough  layer  problem. 

The  experimental  work  was  done  by  acoustic  simulation  tech- 
nique.   A  rough  layer  was  made  from  2  inches  thick  plexiglas 
sheet,  one  surface  of  which  was  machined  to  give  a  layer  with  one 
surface  rough  and  the  other  surface  smooth.   Backscatter  measure- 
ments were  made  for  two  configurations  of  the  layer: 

(1)  the  back  surface  rough  and  the  front  surface  smooth. 

(2)  the  front  surface  rough  and  the  back  surface  smooth. 
The  dependence  of  mean  radar  cross  section  of  layer  <o0> 

(for  both  configurations)  on  the  angle  of  incidence  6  and  o2/X, 
the  ratio  of  the  standard  deviation  of  the  heights  of  the  rough 
surface  to  the  wavelength  of  the  wave  incident  on  the  rough 
surface,  is  plotted  in  Figs.  10  to  13,  Chapter  4.   From  these 
curves  it  may  be  concluded  that  <aQ>   decreases  with: 

(1)  increase  in  roughness  of  layer  (i.e.,  increase  of  oz/X), 

(2)  increase  in  angle  of  incidence. 

This  is  in  agreement  with  the  theoretical  results  of  Chapter  3. 

Various  limitations  of  the  measured  experimental  data  for  the 
rough  plexiglas  layer  by  acoustic  simulation  were :  -''• 

(1)  Plexiglas  is  a  lossy  substance. 

(2)  At  angles  of  incidence  other  than  normal,  the  sound 
energy  transmitted  into  the  plexiglas  layer  propagates 
in  both  the  longitudinal  and  transverse  modes  ;  rather 
than  the  desired  longitudinal  mode  only.   These 
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limitations  are  bound  to  affect  the  measured  results 
while  they  have  been  neglected  in  the  theoretical 
study . 
(3)   The  rough  surface  of  the  layer,  which  was  obtained  by 
machining  it  on  a  shaping  machine,  was  obviously  not 
random  enough,  since  it  was  found  to  have  a  periodic 
component  superimposed  on  a  randomly  rough  component; 
moreover,  the  slopes  of  the  rough  surface  were  not 
distributed  normally  but  had  a  "peak"  at  10° as  well  as  0 
Only  a  partial  attempt  was  made  to  compare  the  experimental 
results  with  the  theoretical  results  and  no  direct  comparison 
could  be  made  due  to  the  limitations  of  the  measured  data,  as 
described  above.   The  experiment  also  violated  the  assumption  of 
the  theory  that  within  the  layer,  the  backscattered  electric 
field  at  the  smooth  front  interface  is  in  the  far  zone  of  dif- 
fraction, i.e., the  mean  depth  of  layer  d»  \   and  d»the  incident 
beamwidth.   The  violation  of  this  assumption  was  unavoidable  be- 
cause: 

(1)  d  could  not  be  increased  beyond  2  inches;  plexiglas 
sheets  thicker  than  2  inches  are  not  easily  available 
commercially.   Even  if  they  were  available,  the  wave 
propagating  within  it  would  be  attenuated  considerably, 
thus  violating  the  assumption  that  the  layer  is  loss- 
less. 

(2)  X  and  beamwidth  could  not  be  decreased  further;  the 
frequency  range  was  limited  due  to  the  limited  number 
of  transducers  and  the  pulsed  oscillator  available. 
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6.   RECOMMENDATIONS  FOR  FUTURE  STUDIES 

In  Chapter  3  the  normalized  radar  cross  section  of  a  rough 
layer  with  the  front  surface  smooth  and  the  back  surface  rough 
has  been  derived  with  many  simplifying  restrictions;  however, 
it  seems  to  be  the  first  such  attempt  to  solve  the  rough  layer 
problem.   The  literature  was  thoroughly  reviewed;  it  seems  that 
no  previous  attempt  has  been  made  to  solve  the  problem  of  back- 
scatter  of  electromagnetic  waves  from  a  rough  layer. 

The  model  to  be  considered  next  should  consist  of  a  rough 
layer  having  the  front  surface  rough  and  the  back  surface  smooth, 
assuming  as  before  that  the  depth  of  layer  is  much  greater  than 
the  incident  wavelength  within  the  layer  so  that  the  back  sur- 
face of  the  layer  may  be  considered  to  lie  in  the  far  zone  of  dif- 
fraction.  It  may  be  pointed  out  that  the  theoretical  analysis 
of  this  model  will  be  much  more  complex  compared  to  the  model 
chosen  in  Chapter  3.   A  statistical  approach  to  this  problem  is 
worth  investigating;  it  might  prove  to  be  less  complex  and  yield 
more  useful  results  as  compared  to  the  physical  optics  method. 
Parks  [30]  gives  a  statistical  solution  for  backscatter  from  sea, 
in  the  near  zone  of  diffraction.   This  approach  may  be  extended 
to  the  case  of  backscatter  from  a  rough  layer.   The  case  of  a 
single  rough  layer  may  be  said  to  be  completely  solved  only  when 
the  solution  is  found  for  a  layer  having  both  the  boundaries 
randomly  rough  and  no  restriction  on  the  mean  depth. 

Meanwhile  an  experimental  approach  may  be  envisaged.   Instead 
of  acoustic  simulation  of  the  electromagnetic  problem,  it  is 
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recommended  that  experiments  should  be  conducted  using  micro- 
waves .   This  will  give  a  greater  flexibility  in  choosing  the 
material  for  the  layer,  and  will  also  yield  more  quantitative 
results.   While  choosing  the  layer  material  for  acoustic  simu- 
lation it  was  observed  that  most  of  the  available  materials 
(solids)  were  lossy  and  had  a  very  high  reflection^  coefficient 
when  the  sound  wave  was  incident  from  water.   Moreover,  the 
energy  transmitted  by  the  longitudinal  acoustic  wave  into  a 
solid  from  a  liquid-solid  interface  has  two  modes  of  propaga- 
tion, the  longitudinal  mode  and  the  transverse  (shear)  mode. 
This  causes  an  error  in  the  observed  results .  At  large  angles 
of  incidence  most  of  the  sound  energy  transmitted  into  the 
solid  propagates  in  the  transverse  mode,  rather  than  the  longi- 
tudinal mode.   Therefore,  care  should  be  taken  not  to  use  large 
angles  of  incidence,  the  limit  depending  upon  the  actual  liquid 
and  solid  used. 

Various  models  of  layer  should  be  made,  their  backscattered 
signal  recorded,  and  their  average  radar  cross  section  calcur 
lated.   The  materials  as  well  as  the  statistics  of  the  target 
may  be  varied.   A  complete  set  of  such  experiments  will  yield 
considerable  information  on  the  statistics  and  electrical  prop- 
erties of  the  layer  material. 
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APPENDIX  I 
Target  Construction,  Properties  &  Statistics 

Target  Construction. 

The  heights  of  the  rough  surface  (one  dimensional)  of  the 
layer  were  chosen  to  belong  to  a  normal  distribution  with 
standard  deviation  oz  =  0.05".   Only  the  area  under  the  normal 
curve  lying  between  -3oz  and  +3oz  was  considered.   These 
heights  were  then  arranged  to  give  a  Gaussian  covariance  func- 
tion (normalized  autocorrelation  function)  with  a  correlation 
distance  of  0.15  inch.   This  particular  correlation  function 
was  obtained  by  trial  and  error  using  a  1620  IBM  computer.   It 
was  noted  that  the  heights  had  to  be  arranged  in  such  a  manner, 
that  there  were  not  many  sudden  variations  in  the  slopes  of  the 
surface  thus  generated.   This  random  curve  of  known  statistics 
was  then  traced  on  the  edge  of  the  plexiglas ,  which  was  then  cut 
by  a  shaping  machine;  care  being  taken  to  follow  the  curve  as 
closely  as  possible. 

After  the  plexiglas  had  been  cut,  the  rough  surface  was 
sampled  and  its  mean  value,  standard  deviation  and  normalized 
autocorrelation  function  was  redetermined.   The  profilometer 
(apparatus  for  determining  the  profile  of  a  rough  surface)  used 
for  sampling  the  target  was  accurate  up  to  one  thousandth  of 
an  inch.   The  sampling  was  done  at  intervals  of  0.05  inch.   The 
measured  standard  deviation  and  correlation  distance  are  0.0465. 
inches  (1.18mm)  and  0.150  inches  (3.81mm),  respectively.   The 
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normalized  autocorrelation  function  is  given  in  Fig.  8,  Chapter 
4.   The  mean  value  was  found  to  lie  almost  midway  between  the 
highest  and  lowest  points  sampled. 

The  estimate  of  the  normalized  correlation  function  was 
calculated  using  the  following  expression: 


N-t 

1 


(hi-m) (hi+T-m) 


i=l 
R(tAx)  ■  


N 
i=l 


where 

hi  is  the  ith  height  in  the  series  of  N  points. 

x   is  the  lag. 

Ax  is  the  sampling  interval. 

N 
_   1  ^S>hi  is  the  sample  mean. 

m   n  <* 

i=l 

At  normal  incidence,  the  diameter  of  the  area  illuminated 
on  the  target  varies  from  1  to  2  inches,  depending  upon  the 
particular  set  of  transducers  used  ,  the  distance  between  the 
transducers  and  the  target  being  32  inches.   Taking  the  mean 
value  of  the  diameter  of  illuminated  area  to  be  1.5  inches,  the 
correlation  distance  of  0.15  inches  was  chosen  to  be  about  one 
tenth  of  this  mean  beam  width.   The  value  of  oz  ■  0.15  inches  was 
chosen  to  give  a  very  rough  surface  for  all  the  four  frequencies 
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used;  the  criterion  for  a  very  rough  surface  being  z>>0.16, 
Tables  II  and  III  list  the  values  of  X  and  ZjL   for  water  and 
plexiglas ,  respectively . 


Table  II.   Wavelength  of  Sound  in  Water. 


Frequency 

in  Mc. 

X  in  mm 

Ratio  oz/X (oz=l 

.18mm) 

1.00 

1.500 

0.787 

1.60 

0.946 

1.250 

2.25 

0.666 

1.770 

3.50 

0.428 

2.760 

The  velocity  of  sound  in  water  was  taken  to  be  1500  meters/sec, 


Table  III.  Wavelength  of  Sound  in  Plexiglas. 


Frequency  in  Mc.  • 

X  in  mm 

Ratio  oz/X (oz=l 

,18mm) 

0.72 

3.860  - 

0.306 

1.00 

2.780 

0.425 

1.28 

2.170 

0.544 

1.60 

1.740 

0.679 

1.90 

1.465 

0.806 

2.25 

1.235 

0.955 

3.00 

0.927 

1.270 

3.50 

0.794 

1.490 

The  longitudinal  velocity  of  sound  in  plexiglas  was  experi- 
mentally determined  to  be  2780  meters/sec. 
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Reflection  Coefficient  at  Water-Plexiglas  Interface;  At  normal 
incidence,  the  acoustic  (characteristic)  impedance  Z  of  a  mater- 
ial is  given  by 

Z  =  pC 

where 

p  ■  density  of  material. 

C  =  longitudinal  velocity  of  sound  in  the 
material. 

For  water,  the  acoustic  impedance  is 

4 
Zw  ■  PlCw  "  15x10 

where 

p^  ■  1.0  gms/cm^ 

Cw  -  15X104  cm/sec. 

For  plexiglas,  the  acoustic  impedance  is 

Zp  ■  p2Cp  =  33.40xl04 

where 

p  2  =  1 • 2  gms/cm 

4 
Cp  =  27.8x10   cm/sec. 

The  reflection  coefficient  r  for  a  plane  sound  wave  incident 
from  water  on  a  water-plexiglas  interface,  is  given  by 

r  ■  ZP  ~  zw 
Zp  +  Zw 

-  0.380  . 
It  may  be  mentioned  that  many  man  hours  were  spent  in  gene- 
rating a  random  curve  of  known  statistics,  constructing  this 
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curve  on  plexiglas,  and  finally  re-calculating  the  statistics 
of  the  rough  surface  thus  obtained.   Two  computer  programs  used 
for  computing  the  mean  value,  variance,  normalized  autocorrela- 
tion function  and  the  radar  cross  section,  respectively,  of  the 
one  dimensionally  rough  surface,  are  given  below. 

COMPUTER  PROGRAM  FOR  DETERMINATION  OF  TARGET  STATISTICS. 

C  MEAN  VARIANCE  AND  AUTOCORRELATION  FUNCTION  OF  ROUGH  SURFACE 

C  FORTRAN  II  LANGUAGE 

C  N  =  TOTAL  SAMPLED* POINTS 

C  N2=  MAXIMUM  SAMPLING  DISTANCE 

C  A (I)  =  SAMPLED  PTS.  OF  ROUGH  SURF ACE, FEED  AT  END  OF  PROGRAM 

99  DIMENSION  K (3000) ,A(3000) 

1  READ  30, N 

30  FORMAT (15) 
READ  31, N2 

31  FORMAT (15) 

2  READ  32,(K(I),I=1,N) 

32  FORMAT (13, 17 14) 
DO  50  1=1 ,N 
A(I)=K(I) 

50  CONTINUE 
AM=N2 
AL=N 
AN=AL-AM 
N1=AN 
STM=0. 
SEM=0. 

3  DO  4  1=1, Nl 
STM=STM+  (A  ( I )  /AN ) 
SEM=SEM+ ( A ( I ) * A ( I ) /AN ) 

4  CONTINUE 
VTR=SEM- (STM*STM) 
DTV=SQRT (ABS (VTR) ) 
PUNCH  12,N 

PUNCH  13,STM 
PUNCH  15, VTR 
PUNCH  14,DTV 
PUNCH  16 
PUNCH  44 

12  FORMAT  (///25HTOTAL  NUMBER  OF  POINTS  IS, 110) 

13  FORMAT (  //10HAVERAGE  =  ,F20.5//) 
15  FORMAT (11HVARIANCE  =  ,F20.5//) 

14  FORMAT ( 2 1HSTANDARD  DEVIATION  =  ,F20.5//) 
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16  FORMAT (///36HNORMALLIZED  AUTO  CORRELATION  FUNCTION) 
44  FORMAT (  //16X,1HT,31X,4HZ (T)/) 

L-0 
25  RES=0. 

DO  8  1=1, Nl 

J=L+I 

RES=RES+ (A ( I ) * A ( J ) ) /AN 

8  CONTINUE 

AUTO- (RES- (STM*STM) ) / (DTV*DTV) 
R=L 

PUNCH  20,R,AUTO 
20  FORMAT(4X,F15.2,6X,F30.10) 

L=L+1 
IF(N2-L)27,25,25 

27  STOP 
END 

> 

COMPUTER  PROGRAM  FOR  CALCULATION  OF  NORMALIZED  RADAR  CROSS 
SECTION  OF  TARGET. 

C      FIND  REFLECTION  COEFFICIENT  AND  NORMALIZED  RADAR  CROSS 

C      SECTION  OF  TARGET  FROM  RECORDED  DATA 

C      FORTRAN  II  LANGUAGE 

C      N  =  SAMPLED  POINTS  FROM  RECORDED  DATA 

C      R  =  REFLECTION  COEFFICIENT 

C      C  -  RECORDED  SIGNAL  FROM  WATER-AIR  INTERFACE 

C      D  -  CHANGE  OF  ATTENUATION  IN  DECIBLES  FOR  RECORDED 

BACKSCATTERED  SIGNAL 
C      SIGMA  =  NORMALIZED  RADAR  CROSS  SECTION 

DIMENSION  A(500) 

DIMENSION  W(10)  ,X(10) 

1  READ  2,N 

2  FORMAT (15) 

READ  3,FREQ,ANGLE 

3  FORMAT (F10. 2, Fl 0.2) 
READ  4,C,D 

4  FORMAT(F10.2,F10.2) 

READ  20, (W(I) ,1-1,5) f(X(I) ,1-1,5) 

20  FORMAT (5A1, IX, 5A1) 
READ  21, (A(I) ,I=1,N) 

21  FORMAT  (F4.0) 
PUNCH  100 

100  FORMAT  (/////)  t    4  ,s 

PUNCH  50,  FREQ, ANGLE, (W(I) ,1=1,5) , (X(I) ,1=1,5) 
50  FORMAT (////12HFREQUENCY  =  ,F7 .2 ,7X,8HANGLE  =  , 
1F7.2,10X,5A1,1X,5A1) 
AN-N 
AVE-0. 
SQAVE-0 . 
FOAVE-0. 
E=LOG(10.)/20. 
D1=D*E 
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A1=C*EXP(D1) 

A2=A1**2 

A4=A2**2 

5  DO  6  1=1, N 
AVE=AVE+(A(I)/AN) 
SQAVE=SQAVE+ (A (I) **2/AN) 
FOAVE=FOAVE+ ( A ( I ) *  *  4/AN ) 

6  CONTINUE 
VARCE=SQAVE- (AVE*  *  2 ) 
DEVIA=SQRT (ABS (VARCE) ) 
VARSQ=FOAVE- (SQAVE**2) 
DEVSQ=SQRT (ABS (VARSQ) ) 
AVE=AVE/A1 
SQAVE=SQAVE/A2 
VARCE=VARCW/A2 
DEVIA=DEVIA/A1 
VARSQ=VARSQ/A4 
DEVSQ=DEVSQ/A2 

PUNCH  16, AVE 

16  FORMAT  (//12HMEAN  OF  R  =  ,  F20.6) 
PUNCH  17,SQAVE 

17  FORMAT (// 16 HMEAN  OF  SIGMA  =  ,   F20.6) 
PUNCH  18, VARCE 

18  FORMAT (// 16 HVARIANCE  OF  R  ■  .F20.6) 
PUNCH  19, DEVI A 

19  FORMAT(//26HSTANDARD  DEVIATION  OF  R  =  ,F20,6) 
PUNCH  2 5, VARSQ 

25  FORMAT (// 20 HVARIANCE  OF  SIGMA  =  ,F20.6) 
PUNCH  26,  DEVSQ 

26  FORMAT (// 3 OHSTANDARD  DEVIATION  OF  SIGMA  =  ,   F20.6///) 
GO  TO  1 

END 
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APPENDIX  II 
Comments  on  Signal  Recording;  Sample  Calculations. 

Comments  on  Recorded  Backscattered  Signal 

Let  the  voltage  recorded  on  the  graphic  level  recorder  be 
V  volts.   Then, 

V  =  klErr 

=  |e  CI) 

r  r 

where,       k  =  a  constant. 

k,  =  the  gain  of  the  electronic  system  between  the 
receiving  transducer  and  the  graphic  level 
recorder, 
r  «  the  distance  of  the  transducers  from  the  target. 
E  ■  the  voltage  proportional  to  the  acoustic  pres- 
sure scattered  by  the  target  (Fig.  14). 
E   =  the  voltage  proportional  to  the  backscattered 

acoustic  pressure  sensed  by  the  receiving  trans- 
ducer (Fig.  14) . 
If  the  reflection  coefficient  of  the  target  is  given  by  R, 

then 

Er  =  REi  C2) 

where,  E.^  is  the  voltage  proportional  to  the  acoustic  pressure 
incident  on  the  target. 
Substituting  (2)  in  (1)  gives 

V  =  -RE.  (3) 

r  * 


83 


v  \\\\v\\  wwwwww.ww 


m  .  u 


o 


ui 
e 


© 


■•• 


Is 

- 
.    .. 


- 

i  i2     q 
x-Jcc 

Q.  LiJ  O 

<  >  o 

C-UJ  Ul 

-.   i- 


03 


CO 


L  . 

- 

:: 


^ 


CO 

-1 
o 

:- 
< 
_i 
z> 
o 
_J 
< 
o 

X 
hi 

H 

< 
O 
CO 

o 

< 
_ 


_ 


84 


when  R  =   1,    (3)    reduces  to 


V,    =  £e.  (4) 

1       r  i 

Dividing  (3)  by  (4)  gives 

»-*!     ■  «) 

It  is  well  known  that 

R2  =  aQ  (6) 

Substituting  (6)  in  (5)  gives 
V2 

Therefore,  if  the  square  of  the  backscattered  signal  record- 
ed by  the  graphic  level  recorder  for  any  target  is  divided  by  the 
square  of  the  backscattered  signal  from  a  water-air  interface 
(for  which  R  =  1)  at  the  same  distance,  it  will  give  the  normal- 
ized radar  cross  section  of  the  target. 

Sample  Calculations 

Target:   Layer  with  a  rough  front. 

Backscattered  signal:   From  rough  front  only. 

f  -  1.0  Mc. 

e  =  0° 

V.=  0.66  volts,  (water-air  interface) 

Total  attenuation  present  (for  V]_  =  0.6  volts)  in  the  receiving 

electronic  system  =  65  db. 

Root  mean  square  of  backscattered  signal  from  layer  /(<V^>)  ■  0.286 

volts 
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Total  attenuation  present  (for  /(<V2>)  =  0.286  volts)  in  the 
receiving  electronic  system  =  50  db. 

Therefore,  the  increase  in  the  gain  of  electronics  system  to 
record  the  backscattered  signal  from  layer  =  65-50 

=  15  db. 
Since  the  gain  of  the  receiving  electronics  system  was  increased 
while  recording  the  signal  backscattered  from  the  layer,  the 
value  of  V1  will  have  to  be  modified.   Denote  this  modified  value 
by  V^   .   Then, 


20  log  Xi  =  15  db 


or 


vl     15 
lo9  0766  "  20  "  °'75 


or 


vy  =  0.66x5.62 
=  3.72  volts  . 
The  value  of  <oQ>  is  given  by 
f2^     fn_2flfi}2 


<V2>     (0.286)*    cq     -4 
o>  "  TVTTZ"  =  (3:72)2   =  59.2x10 


<o«>  = 
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APPENDIX  III 
Tables  of  Measured  Data 


Table  IV.   Measured  data 

Cor  rough  back  configuration  of  layer. 

Frequency 
in  Mc, 

Ratio  oz/X 

Angle  of 
incidence 
in  degrees 

<a0> 

Standard 
deviation  of 
<o0> 

0.72 

0.306 

0 

5 

10 

0.008632 
0.008300 
0.008893 

0.006411 
0.005937 
0.005540 

1.00 

0.425 

0 

5 

10 

20 

0.002217 
0.001208 
0.001666 
0.000340 

0.002229 
0.001078 
0.001995 
0.000288 

1.28 

0.544 

0 

5 

10 

0.001203 
0.001116 
0.001249 

0.000660 
0.000867 
0.001040 

1.60 

0.679 

0 

5 

10 

20 

0.000480 
0.000638 
0.001018 
0.000221 

0.000427 
0.000674 
0.000957 
0.000190 

1.90 

0.806 

0 

5 

10 

0.000166 
0.000168 
0.000269 

0.000114 
0.000131 
0.000241 

2.25 

0.955 

0 

5 

10 

20 

0.000111 
0.000158 
0.000273 
0.000046 

0.000116 
0.000132 
0.000210 
0.000046 

3.00 

1.270 

0 

5 

10 

0.000093  * 

0.000014 

0.000019 

0.000044 
0.000008 
0.000010 

3.50 

1.490 

0 

5 

10 

20 

0.000054 

0.000009 

0.000004 

1  0.000000 

0.000079 
0.000007 
0.000004 
0.000000 
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Table  V.   Measured  data 

for  rough  front  configuration  of  layer 

Frequency 
in  Mc. 

Ratio  az/\ 

Angle  of 
incidence 
in  degrees 

<o0> 

Standard 
deviation 
of  <o0> 

1.00 

0.787 

0 

5 

10 

20 

0.006088 
0.003782 
0.003544 
0.001468 

0.006694 
0.006298 
0.006640 
0.002024 

1.60 

1.250 

0 

5 

10 

20 

0.001422 
0.000597 
0.001910 
0.000423 

0.001497 
0.000454 
0.002149 
0.000590 

2.25 

1.770 

0 

5 

10 

20 

0.001353 
0.002490 
0.002564 
0.001067 

0.000852 
0.001563 
0.002345 
0.001899 

3.50 

2.760 

0 

5 

10 

20 

0.000651 
0.000680 
0.000715 
0.000388 

0.000424 
0.000417 
0.000487 
0.000589 
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APPENDIX  IV 
Scattering  from  Periodically  Rough  Surfaces. 

The  equation  from  grating  theory,  given  below,  indicates 
that  a  periodic  surface  scatters  the  incident  wave  into  a  num- 
ber of  different  modes  (or  lobes)  as  compared  to  a  rough  surface 
which  scatters  diffusely. 

sine,.  =  sinen  +  m  X     (m  =  0,  +  1,  +  2,...)      (1) 
2m  A 

where 

6,  ■  angle  of  incidence 

X  =  wavelength  of  incident  wave-* 

A  =  wavelength  of  periodic  surface 

m  =  any  integer 

6_  =  scattering  angle  corresponding  to  each  mode. 

The  scattering  angles  62m  are  determined  by  the  grating 
equation  (1) .   To  each  integer  m  there  corresponds  a  scattered 
mode  propagated  in  the  direction  82m«   The  total  number  of  pos- 
sible modes  is  limited  by  the  condition. 

|sin02m|<  1  (2) 

The  mode  m  =  0  is  seen  from  (1)  to  be  the  specular  mode. 
The  modes  m  =  +  1  lie  to  either  side  of  the  specular  direction. 
The  modes  continue  thus  to  either  side  of  the  specular  mode  until 
the  last  modes  that  will  satisfy  (2)  are  reached. 

If  X/a  is  small,  it  follows  from  (1)  that  m  will  run  through 
a  large  number  of  integral  values  before  (2)  is  violated,  so  that 
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if  the  wavelength  of  the  incident  radiation  is  small  compared 
to  the  period  or  wavelength  of  the  surface,  the  incident  wave 
will  be  broken  up  into  many  scattered  waves. 
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ABSTRACT 

The  radar  cross  section  of  a  rough  surface  is  often  used  to 
predict  the  dielectric  and  statistical  properties  of  the  scatterer, 
In  an  attempt  to  extend  this  approach,  a  target  consisting  of  a 
rough  layer  with  the  front  surface  smooth  and  the  back  surface 
statistically  rough  was  chosen  for  this  study.   Theoretical  ex- 
pressions were  derived  for  the  mean  value  of  the  radar  cross  sec- 
tion <o0>  of  the  rough  layer.   These  expressions  for  <o0>  turn 
out  to  be  complicated  functions  of  the  angles  of  incidence  and 
scatter,  frequency,  mean  depth  of  layer,  dielectric  and  statis- 
tical properties  of  the  rough  layer,  etc. 

To  better  understand  the  behaviour  of  a  rough  layer,  exper- 
imental work  was  conducted  using  the  acoustic  simulation  tech- 
nique.  A  plexiglas  layer  48  inches  long,  6  inches  wide  and  2 
inches  thick,  was  carefully  machined  on  one  surface  to  obtain  a 
statistical  variation  in  the  thickness  with  nearly  Gaussian  dis- 
tribution and  nearly  Gaussian  correlation  function.   The  layer 
was  placed  in  a  water  tank  with  water  on  one  side  and  air  on  the 
other,  and  was  struck  by  pressure  waves  from  electro-acoustic 
transducers  of  different  frequencies,  in  the  usual  manner. 

The  backscattered  signal  was  recorded  for  eight  frequencies, 
four  angles  of  incidence  and  two  different  configurations  of  the 

layer: 

(a)  The  front  surface  smooth  and  the  back  surface  rough. 

(b)  The  front  surface  rough  and  the  back  surface  smooth. 
From  the  measured  data,  the  following  curves  were  drawn: 


(1)  <o0>  versus  oz/\,    the  back  interface  being  rough. 

(2)  <cjq>  versus  angle  of  incidence,  the  back  interface 
being  rough. 

(3)  Same  as  (1),  but  with  the  front  interface  rough. 

(4)  Same  as  (2) ,  but  with  the  front  interface  rough. 

Here  o,  and  X  are  the  standard  deviation  of  the  heights  of  the 
z 

rough  surface  of  layer,  and  the  wavelength  of  the  incident  wave, 
respectively. 

Only  a  partial  attempt  was  made  to  check  the  experimental 
results  with  those  obtained  from  theory,  due  to  certain  limita- 
tions of  the  experimental  data.   The  theoretical  as  well  as 
experimental  results  indicate  that  the  average  radar  cross  sec- 
tion <oQ>  of  a  layered  target  decreases  as  the  rough  surface  of 
layer  becomes  more  rough  (i.e.  az/\   increases) ,  and  as  the  angle 
of  incidence  becomes  more  oblique.   Other  implications  of  the 
theoretical  expressions  and  the  measured  curves  are  also  dis- 
cussed. 


